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ABSTRACT 


The  state  of  the  art  In  steady-state  augmentors  and  jet  pomps  la  briefly  re¬ 
viewed  and  a  general  performance  theory  developed.  Generalized  charts  are 
presented  giving  the  augmentation  ratio  obtainable  from  an  optimized  eductor, 
together  with  the  associated  geometrical  and  fluid-dynamic  parameters.  This 
theory  is  shown  to  give  good  agreement  with  experiment.  Experimental  mea¬ 
surements  made  by  various  investigat  )rs  in  the  past  do  not  achieve  the  pre¬ 
dicted  optimum  performance,  however,  because  of  various  deviations  from 
optimum  design  in  their  test  eductors. 

The  performance  of  eductors  in  an  axial  stream  and  the  total  head  rise  obtain¬ 
able  through  an  eductor  are  both  Investigated  theoretically.  The  report  con¬ 
cludes  with  a  brief  survey  of  various  unconventional  eductors,  including  multi¬ 
stage  units.  Coanda  eductors  and  crypto-steady  flow  devices. 
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FOREWORD 


The  work  herein  reported  was  carried  out  by  Peter  R.  Payne,  Inc. ,  Rockville, 
Maryland,  in  compliance  with  U.  S.  Army  Aviation  Materiel  Laboratories  Contract 
No.  DA  44-177-AMC-337(T). 

The  principal  investigator  and  report  author  was  Mr.  Peter  R.  Payne.  Impor¬ 
tant  contributions  were  made  by  Messrs.  James  O.  Justice  and  Alastair  Anthony. 

In  addition  to  work  carried  out  under  the  above  contract,  this  report  contains,  in 
Chapter  Nine,  the  results  of  investigations  made  for  the  Office  of  Naval  Research 
under  Contract  No.  Nonr-4626<00).  Because  of  its  appropriateness  to  the  present 
report,  Mr.  Ralph  D.  Cooper,  Head,  Fluid  Dynamics  Branch,  Office  of  Naval  Re¬ 
search,  approved  the  inclusion  of  this  work  in  this  publication.  The  report  also 
contains  material  generated  in  the  course  of  a  company  funded  research  program, 
much  of  which  was  originally  reported  in  AIAA  Paper  No.  64-798,  which  appears 
as  Reference  28. 
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INTRODUCTION 


This  re,  rt  is  concerned  with  reviewing  the  state  of  the  art  in  jet  augmentation. 
An  augmentor  is  broadly  defined  as  a  device  which  in  some  way  augments  a  jet  - 
usually  referred  to  as  the  primary  jet  -  by  allowing  ambient  fluid  to  mix  with  it 
to  achieve  a  desired  result.  Since  there  appears  to  be  no  universally  accepted 
terminology,  the  following  definitions  are  used  in  this  report. 

Eductor  (or .occasionally.  Inductor) 

A  generic  term  describing  devices  in  which  a  primary  fluid  or  gaseous  jet  gives 
up  energy  to  a  secondary  flow.  The  exit  flow  thus  has  a  lower  velocity  than  the 
primary  jet,  but  the  mass  flow  is  increased.  The  primary  and  secondary  fluids 
are  not  necessarily  the  same. 

Injector  (for  which  an  alternative  is  Jet  Pump) 

An  injector  is  an  eductor  used  to  pump  fluid  against  a  back-pressure,  such  as 
feed  water  to  a  steam  boiler,  for  example,  or  air  into  the  plenum  chamber  of 
a  GEM.  An  injector  generally  has  a  larger  diffuser  area  ratio  than  the  equiva¬ 
lent  thrust  augmentor  of  the  same  mass  flow  ratio. 

Augmentor 

A  thrust  augmentor  is  an  eductor  which  is  optimized  for  maximum  thrust  in¬ 
crease;  a  mass  flow  augmentor  is  designed  for  maximum  mass  flow  increase. 

In  practice  there  is  little  difference  between  the  two. 

Ejector 

An  ejector  may  be  defined  as  an  eductor  which  exhausts  air  to  ambient  from  a 
region  of  lower  pressure,  as  in  ’’pumping  out"  a  vacuum  chamber,  for  example. 
Although  closely  analogous  to  an  injector,  it  has  significant  geometrical  differ¬ 
ences  if  the  mixing  section  pressure  is  very  low. 

The  basic  geometry  common  to  ail  such  devices  is  shown  in  Figure  1. 
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In  its  most  commonly  used  application,  an  eductor  augments  the  mass  flow  a 
jet  in  order  to  increase  its  thrust.  In  this  application  it  functions  as  a  "trans¬ 
former,  "  in  exactly  the  same  way  as  a  pressure  jet  helicopter  rotor,  for  ex¬ 
ample,  transforms  the  high  energy,  low  mass  flow  gas  energy  from  a  turbine 
exhaust  to  a  high  mass  flow,  low  energy  flow  through  the  rotor. 

The  eductor  stands  alone  among  other  gas  energy  transformers  in  that  it  re¬ 
quires  no  moving  parts  to  "transform"  high  pressure,  low  mass  flow  gas  energy 
to  a  low  pressure,  high  mass  flow.  This  advantage  is  potentially  of  great  impor¬ 
tance.  Unfortunately,  as  shown  in  Figure  2,  the  "transformer  efficiency"  of  ex¬ 
isting  eductors  is  very  low  indeed,  and  it  more  than  offsets  their  savings  in 
weight  and  mechanical  complexity.  But,  if  by  properly  conducted  research  the 
ability  to  reliably  predict  eductor  performance  can  be  achieved,  it  may  be 
possible  to  discover  ways  of  upgrading  the  overall  efficiency  to  the  point  where 
eductors  are  competitive  for  many  applications. 
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The  Efficiency  of  a  Static  Eductor  in  Comparison  With  Some 
Alternative  Methods  of  Transforming  Gas  Energy  into  Static  T'  rust. 


Chapter  1 


A  BRIEF  REVIEW  OF  THE  LITERATURE 


The  eductor  principle  has  been  employed  perhaps  longer  than  can  be  traced, 
and  has  certainly  occupied  the  attention  of  many  modern  engineers  and  inven¬ 
tors.  Timothy  Hackworth  and  George  Stevenson  employed  eductors  to  obtain  a 
forced  draft,  giving  more  efficient  combustion  in  early  railway  locomotives, 
circa  1830.  This  led  naturally  to  steam-powered  water  jet  pumps  for  the  pur¬ 
poses  re-charging  boilers  under  pressure  and  to  the  fairly  extensive  em¬ 
ployment  of  the  eductor  principle  through  the  field  of  mechanical  engineering. 

Naturally  enough, these  early  developments  were  based  on  "cut  and  try"  tech¬ 
niques.  Efficiencies  were  very  low.  Even  today,  "Marks  Handbook, "  a  uni¬ 
versally  recognized  reference  for  the  mechanical  engineer,  says  ", . .  the  ex¬ 
ternal  work  is  usually  about  two  percent  of  the  heat  given  up  by  the  steam. " 

Early  attempts  at  analysis  were  based  upon  the  energy  equation.  Flugel  (Ref¬ 
erence  1)  summarizes  much  of  this  early  work,  in  addition  to  making  a  consid¬ 
erable  number  of  original  and  valuable  contributions.  In  particular  he  identified 
the  fact  that  there  is  an  "optimum"  mixing  pressure,  and  drew  attention  to  the 
fact  that  the  mixing  section  cannot  be  accurately  designed  by  means  of  one¬ 
dimensional  flow  theory.  To  quote  Reference  1: 

"The  greatest  drawback  of  the  energy  equation  is  that  the  assumption  of  complete 
intermingling  of  driving  fluid  and  delivered  fluid  at  constant  pressure  Pj  is  ut¬ 
terly  incorrect,  resulting  from  the  fact  that  the  narrowest  section  of  the  mixing 
nozzle  must  be  designed  from  the  'theoretically*  narrowest  section.  In  conse¬ 
quence,  the  calculation  method  does  not  give  the  correct  information  about  the 
most  important  section  of  a  jet  pump  and  no  data  whatsoever  regarding  the  nec¬ 
essary  leRgth  ratios. 

"Practice  proves  that,  under  normal  operating  conditions,  the  intermingling  is 
practically  always  accompanied  by  a  pressure  rise*,  so  that  at  the  end  of  the 
mixing  nozzle  usually  a  substantial  or  even  a  major  part  of  the  pressure  rise 
will  already  have  been  achieved.  This  fact  gives  consideration  to  the  subse¬ 
quent  method  largely  evolved  from  the  impulse  theorem  according  to  which, 
first  of  all,  a  length  for  the  mixing  nozzle  is  assumed  necessary  to  assure  ade¬ 
quate  mixing.  The  assumption  of  constant  pressure  Pj  during  mixing  can  be 


*  That  is,  when  the  mixing  takes  place  in  a  constant  area  duct. 
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dropped  again  later.  The  mixing  advances  in  the  narrowest  section  of  the  mix¬ 
ing  nozzle  under  a  pressure  rise  until  a  practically  perfect  intermingling  has 
been  achieved  at  the  end  of  the  narrowest  section.  This  part  of  the  pressure 
rise  is  readily  computed  according  to  impulse  theorem.  It  is  clear  and  closely 
according  to  practice  that  with  the  comparative  slowness  of  the  mixing  process 
the  length  L  of  the  narrowest  section  obviously  must  be  fairly  great  in  order 
to  assure  adequate  mixing  (which,  in  turn,  is  the  premise  for  achieving  proper 
energy  conversion  in  the  adjoining  diffuser);  it  is  therefore  recommended  that 
L  -  10  d  .  It  may  be  stated  that  very  favorable  results  have  already  been 
achieved  with  jet  pumps  computed  on  this  basis,  which  are  also  in  good  agree¬ 
ment  with  the  preliminary  calculation,  as  will  be  reported  elsewhere.  " 

Flugel's  work  is  apparently  not  toe  well  known,  despite  the  fact  that  it  is  a  very 
exhaustive  study,  and  illuminates  a  number  of  problem  areas  which  have  troubled 
subsequent  investigators 

The  first  American  workers  to  develop  a  theoretical  treatment  of  eductors  were 
McClintock  and  Hood  (Reference  3).  Starting  from  essentially  the  same  assump¬ 
tions  as  the  theory  of  the  present  report,  they  obtained  an  equation  for  thrust 
augmentation  ratio  which  is 

J2/J0  =  (1  +  n)2/(A2/Aj)  (1) 


for  one-dimensional  flow,  with  equal  primary  and  secondary  densities.  This  is 
merely  another  way  of  writing  Newton’s  second  law,  of  ccrse,  since 


total  momentum  flux  out  of  eductor  •  ,  ..  • 

^2^0  ~  momentum  flux  of  primary  jet  “  u2mj'  +  n  u<imj 

and  un  =  mn/cAn  • 

Thus,  in  the  McClintock  and  Hood  analysis,  a  solution  must  be  found  by  obtain¬ 
ing  a  relationship  between  the  entrainment  ratio  (  n  )  and  the  geometrical  prop¬ 
erties  of  the  eductor.  This  was  done  by  introducing  an  empirical  constant  which 
was  determined  from  eductor  tests.  Additionally, the  augmentation  equation  was 
multiplied  by  an  empirical  "thrust  coefficient,  ”  nominally  intended  to  allow  for 
deviations  from  one-dimensional  flow.  The  disadvantage  of  introducing  some 
empiricisms  into  a  theory  is  obvious.  It  is  as  well  to  note  that  McClintock 
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and  Hood  were  concerned  primarily  with  the  use  of  eductors  to  promote  piston 
engine  cooling. 

Shortly  after  the  advent  of  the  jet  engine, several  theoretical  and  experimental 
studies  of  very  low  entrainment  ratio  ’’ejectors"  were  made  (References  4  and 
5).  This  stemmed  from  the  need  to  cool  the  inside  of  jet  engine  nacelles,  Ob¬ 
viously  the  jet  engine  exhaust  provided  a  convenient  means  of  doing  so.  The 
results  of  these  investigations  are  of  littie  value  in  the  broader  context  of  high 
entrainment  ratio  eductors  intended  to  give  augmentation,  however. 

In  1949,  von  Karman  (Reference  6)  treated  the  case  of  a  loss-free  eductor  with 
a  parallel  wall.  For  area  ratios  A/Aj>15  and  one-dimensional  flow, his  results 
are  given  by  the  approximate  equations 


uju  .  *(A  ./A)*  -  A  ./A  (2) 

^  J  J  J 

JjJ  =  (1  +  A/A,)  UkJA)1*  -  A  ./A)2  .  (3) 

*  °  V  J  J 

Von  Karman  pointed  out  that  this  implied Jp/Jo<2.0  t  always,  so  long  as  the  flow 
was  one-dimensional.  This  of  course  applies  only  to  constant  diameter  eductors, 
higher  values  being  obtainable  when  a  diffuser  follows  the  mixing  section. 

He  also  pointed  out  that  non-uniform  inlet  flow  could  result  in  greater  augmenta¬ 
tion,  an  observation  already  implied  by  the  work  of  McClintock  and  Hood.  This 
important  observation  does  not  seem  to  have  been  pursued  by  subsequent 
workers. 

Von  Karman’s  analysis  applies  only  to  rectilinear  flow,  because  Bernoulli’s 
equaticn  is  used  to  relate  local  speed  to  local  static  pressure.  It  is  not  easy 
to  see  low  a  non-uniform  rectilinear  flow  could  be  developed  without  total  pres¬ 
sure  loss.  This  may  be  the  reason  for  the  apparent  neglect  of  his  suggestion. 
Never;heless,  a  more  thorough  analysis  of  this  aspect  appears  well  worthwhile. 

It  is  of  interest  to  note  that  augmentation  ratios  greater  than  the  theoretical 
value  of  2.0  were  achieved  experimentally  by  Payne  (Reference  11)  in  1954, 
using  constant  diameter  ducts. 

A  number  of  the  larger  aerospace  companies  have  conducted  experimental  stud¬ 
ies  since  the  publication  of  these  and  other  early  theoretical  studies.  Generally 
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speaking,  they  made  no  effort  to  extend,  or  even  consolidate,  the  theory,  but 
concentrated  on  experimental  measurements.  As  might  be  expected,  the  re¬ 
sults  obtained  with  the  relatively  large  and  sophisticated  test  rigs  were  no 
better  than  the  results  obtained  in  small-scale  laboratory  tests  by  other  ex¬ 
perimenters.  From  hindsight  it  is  easy  to  identify  where  deviations  from 
optimum  design  occurred  in  these  large-scale  programs.  In  the  absence  of 
adequate  theoretical  work,  however,  it  was  obviously  difficult  for  the  in¬ 
vestigators  to  avoid  these  deviations  at  the  time  they  were  conducting  their 
programs. 

In  1955,  Szczeniowski  (Reference  12)  published  a  one-dimensiona!  flow  anal¬ 
ysis  of  the  eductor.  He  considered  a  quite  general  case  of  varying  mixing 
chamber  geometry,  using  the  solutions  for  constant  pressure  and  constant 
area  mixing  as  particular  solutions  of  the  general  case.  He  discussed  the 
mathematical  problem  of  determining  the  optimum  static  pressure  distribution 
along  the  mixing  chamber  without  reaching  any  firm  conclusions.  No  absolute 
optimum  exists,  but  some  pressure  distributions  are  better  than  others. 

Szczeniowski's  analysis  did  not  include  the  effects  of  internal  duct  losses,  which 
subsequent  research  has  shown  to  be  of  dominating  importance. 

An  important  program  of  research  was  carried  out  by  Helmbold  (References  9 
and  13)  during  1953  and  1954,  in  which  he  studied  various  aspects  of  the  eductor 
cycle.  Of  particular  interest  was  his  design  of  a  constant  static  pressure 
mixing  section,  which  gave  31  percent  greater  overall  efficiency  than  an  equiv¬ 
alent  cylindrical  mixing  chamber.  This  is  in  good  agreement  with  Helmbold's 
theoretical  predictions.  It  would  appear  that  the  design  and  incorporation  of  a 
constant  static  pressure  mixing  section  in  the  XV-4A  research  aircraft  eductor 
would  have  provided  worthwhile  gains  in  the  Bystem. 

Recently,  Hill  (Reference  56)  effected  a  major  consolidation  of  the  theoretical 
aspects  of  the  mixing  chamber  phenomena.  This  work  may  be  represented  as 
a  major  step  toward  the  goal  of  a  consistently  reliable  predictive  ability. 

Wells  (Reference  15)  advanced  Helmbold's  constant  mixing  pressure  c  anber 
work  a  step  further  by  modifying  it  to  maintain  a  favorable  static  pressure 
gradient  as  far  as  possible,  and  hence  to  reduce  skin  friction  by  increasing  the 
length  of  the  laminar  boundary  layer.  He  also  applied  boundary  layer  control, 
in  the  form  of  blowing  from  a  slot,  at  the  diffuser  throat.  Although  exact  com¬ 
parison  with  Helmbold's  experiment  is  difficult,  there  seems  little  doubt  that 
these  innovations  resulted  in  a  further  improvement  in  efficiency. 
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Although  this  brief  review  of  the  literature  is  by  no  means  exhaustive,  it  seems 
clear  that  a  great  deal  of  research  work  remains  to  be  done  before  either  ade¬ 
quate  predictive  ability  is  achieved  or  maximum  theoretical  efficiency  can  be 
approached  in  new  designs.  It  seems  equally  clear  that  hardware  applications 
of  the  eductor  principle  have  not  profited  to  the  extent  made  possible  by  existing 
research  work.  Further  significant  improvements  can  be  made,  provided  ex¬ 
periment  and  theory  are  allowed  to  develop  side  by  side  in  the  laboratory. 
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Chapter  2 


SIMPLE  PERFORMANCE  RELATIONSHIPS  FOR  A  PROPULSOR 


A 


2 


Figure  3.  A  Generalized  Propulsor. 


IDEAL  PROPULSIVE  EFFICIENCY  DEPENDS  ONLY  ON  THE  JET  VELOCITY 

Consider  the  generalized  propulsor  in  Figure  3,  which  accelerates  ambient  fluid 
from  a  free -stream  velocity  uQ  to  a  final  velocity  of  .  The  temperature 
of  the  downstream  jet  is  not  necessarily  ambient,  so  that  the  efflux  density  p^p^. 

The  exit  mass  flow  is  (  p^2u2  )*  so  that  its  kinetic  energy  power  is 


H(pzA2U2)  u 


(4) 


Since  the  fluid  originally  had  some  kinetic  energy  before  entering  the  propulsor, 
the  amount  of  energy  added  by  the  propulsor  is 


APJ  =Jf(f>2A2u2)  (u2  -  uo2)  * 


(5) 
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The  developed  thrust  (  T  )  is  equal  to  the  rate  of  change  of  momentum.  Thus, 
in  unaccelerated  motion, 


T=  (?2A2u2)  (u2  ~  uo>  •  (6) 


The  ideal  efficiency  is  defined  as  the  ratio  of  useful  propulsive  power  to  the  power 
lost  in  the  jet.  That  is, 

\  =  KV APj  =  q(u2  -  uQ)  /  (v\  -  up 


*  2/  1  +(u2/u0}  ’  (?) 

Thus ,  the  ideal  efficiency  is  a  function  only  of  the  velocity  ratio  u  p/u ,  and  has 
a  maximum  of  unity  when  up  -  uQ  t  the  zero  thrust  condition.  The  temperature 
of  the  propulsor  efflux,  which  would  cause  a  change  in  density  of  the  fluid,  has  no 
influence  upon  the  efficiency. 

THE  RELATIONSHIP  BETWEEN  THRUST  AND  TOTAL  HEAD  INCREASE  IN 
THE  WORKING  FLUID 

From  Figure  3, the  total  head  of  the  working  fluid  is  increased  by 
A#  =  hP2'u2  ~  ouo 

as  it  passes  through  the  propulsor.  Using  this  relationship  to  define  Up  in 
terms  of  u0  and  AH  >  the  thrust  equation  becomes 

T/qoU1A2  =  (1  ~  u0/u2)  *  (8) 


Writing  AH  =  &H/%pu%  also  gives 


(9) 


LH  =  fi2/o0  (u2/u0)2  -  2  ; 

therefore, 

T/onu2A  2  ~  (1  +  tA)  -  (po/p0^  d  +  • 


For  the  static  case  {  uQ  -  0  ),this  reduces  to 


T  =  2UJ  A2 


(10) 


That  is  to  say,  changes  in  density  of  the  working  fluid  do  not  influence  the  result, 
but  only  the  total  head  rise.  For  a  propeller,  where  A2  ~  under  static 


conditions  f 


T  =  LB.  it/?2  , 


(11) 


a  result  which  is  already  familiar.  The  ramjet, on  the  other  hand  .is  a  device 
which  cannot  increase  total  head,  but  must  obtain  thrust  by  reducing  the  density 
of  the  efflux  by  addition  of  heat  in  the  case  of  an  air  ramjet,  and  by  introducing 
gas  bubbles  into  a  water  ramjet.  In  the  ideal  (loss  free)  case,  where  LH  =  0  t 
the  thrust  is  then  given  by 


fpg/p^  . 


(12) 


The  eductor  normally  operates  under  static  (  uo  ~  0  )  conditions,  so  that  at¬ 
tention  can  be  focused  on  its  ability  to  raise  the  total  head  of  the  entrained  fluid. 
The  temperature  of  the  primary  jet  is  therefore  unimportant. 

When  the  eductor  is  moving  through  a  fluid,  however,  a  hot  primary  jet  could 
result  in  an  advantageous  increase  in  efficiency,  since  it  would  result  in  p2<^o  . 


Chapter  3 


VISCOUS  MIXING 

Viscosity  is  a  measure  of  a  fluid's  capacity  to  carry  a  shear  stress.  When  the 
stress  is  low,  most  fluids  react  in  a  "Newtonian"  manner  .where  shear  stress 
varies  linearly  with  strain,  the  relationship  being 


T  =  U  3 u/llj  , 


(13) 


or  (stress)  =  (coefficient  of  viscosity)  x  (velocity  gradient). 

When  the  velocity  gradient  exceeds  a  certain  critical  value  usually  defined  in 
terms  of  a  "critical  Reynolds  number,"  the  fluid  no  longer  supports  shear 
stress  in  a  laminar  manner.  Instead,the  smooth  flow  breaks  up  into  eddies 
and  vortices,  giving  what  is  loosely  referred  to  as  turbulent  flow,  tt  seems 
likely  that  the  size  of  these  apparently  random  eddies  decreases  with  increas¬ 
ing  Reynolds  number ,  but  the  whole  picture  is  still  very  poorly  understood. 

The  concept  of  an  "effective  viscosity  coefficient"  is  still  useful  when  the  flow  is 
turbulent.  Its  value  is  much  greater  than  the  laminar  value  (  u  )  for  the  same 
fluid,  and  currently  there  are  only  empirical  expressions  for  the  effective  turbu¬ 
lent  flow  viscosity,  which  differ  widely  in  their  macroscopic  detail.  The 
earliest  of  these  was  Prandtl’s  "mixing  length"  hypothesis, 


t  =  pl2\*u/?y\ (*u/dy)  , 


(14) 


*  being  the  "mixing  length"  and  p  the  fluid  mass  density.  In  the  free  tur¬ 
bulent  shear  flow  characteristic  of  a  jet  or  wake,  it  is  usual  to  assume 


l  *  , 


(15) 


where  c  is  an  empirical  constant  that  is  obtained  from  experimental  data,  and 
the  exponent  (  m  )  is  evaluated  from  the  appropriate  requirements  (usually  "sim¬ 
ilarity"  requirements)  of  the  particular  problem.  (Taylor's  "vorticity  transport 
model,"  of  more  recent  date,  also  leads  to  Equation  14.) 
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Equation  (14)  can  be  criticized  on  a  number  of  grounds,  and  later  workers  have 
developed  alternative  formulations;  notably  the  following: 


The  'Virtual  Kinetic  Viscosity  £  11 


t  =  pe  3u/3y  f e  =  const.)  . 


(16) 


The  "Constant  Shear  Coefficient"  Model  of  Prandtl 


t  =  pe  'du/'dy 


c  =  kb  (u  ~u  .  ) 
max  vn-n  , 


where  k  -  an  empirical  constant  obtained  from  experiment,  and 
width  of  the  mixing  zone. 


(17) 


^  is  the 


Vot  Karman’s  "Similarity"  Hypothesis  gives 


t  =  pk2  ( du/ty)  ^  |  du/'dy  \/ (32u/3j/2)  2 


(18) 


and  finally,  the  "interm:  ttency  model"  of  Townsend,  in  which  Equation  (16)  is 
multiplied  by  an  intermittency  factor 


y  (y>  du 


(19) 


(a  )  and  (  s  )  being  empirical  constants. 

Despite  this  wide  range  of  hypotheses  (  and  there  are  several  more),  all 
give  about  the  same  v  jIoc  ity  distribution  when  used  to  solve  viscous  mixing 
problems. 

•«» 

SOME  PRACTICAL  ASPECTS  OF  VISCOSITY 

\ 

Because  of  viscosity,  a  true  velocity  discontinuity  can  never  exist  between  two 
streams  of  fluid.  As  soon  as  an  air  jet  emerges  from  a  nozzle,  for  example,  it 
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starts  to  carry  along  with  it  some  of  the  previously  stationary  atmospheric  air, 
and  in  the  process  of  accelerating  this  "entrained"  air,  gives  up  some  of  its  own 
momentum.  For  this  reason, the  velocity  profile  of  a  jet  changes  rapidly  as  it 
travels  downstream  from  its  nozzle,  and  accelerates  progressively  more  and 
more  of  the  free  air  around  it,  until  it  soon  loses  ail  resemblance  to  a  discrete 
jet. 

Also,  since  the  ~ree  air  in  the  immediate  vicinity  is  accelerating  to  move  with 
the  jet,  the  pressure  distribution  over  any  solid  bodies  near  the  jet  will  be  mod¬ 
ified  by  the  secondary  airflow,  the  results  of  r.-hich  may  be  beneficial  or  adverse, 
depending  on  the  nature  of  the  problem. 

The  best  known  example  of  this  effect  is  the  "mixing  drag"  experienced  by  a  jet 
engine  in  an  aircraft.  Since  the  free  air  is  accelerated  by  the  jet,  an  area  of 
suction  is  generated  over  the  rear  portions  of  the  nacelle.  The  forces  caused 
by  this  suction  are  inclined  backwards  because  of  the  shape  of  the  nacelle  and  their 
horizontal  component  constitutes  a  drag  which  reduces  the  total  thrust  of  the 
engine -nacelle  combination.  In  this  case  the  mixing  drag  forces  are  usually 
quite  small,  amounting  to  perhaps  1  percent  or  2  percent  of  the  engine  thrust. 
However,  it  is  easy  to  see  that  in  the  case  of  a  long  two-dimensional  jet  whose 
periphery  might  be  ten  times  the  periphery  of  the  equivalent  circular  jet,  the 
jet  drag  could  amount  to  10  -  20  percent  of  the  total  jet  momentum  flux 
under  the  same  geometrical  and  flow  conditions.  Clearly,  this  would  have  a 
large  effect  upon  the  determination  of  optimum  jet  thickness,  and  a  significant 
effect  upon  the  calculated  performance  of  an  aircraft  which  used  a  jet  flap  for 
propulsion,  or  on  the  performance  of  an  annular  jet  GEM. 

Kuchemann  (Reference  83)  hat  observed  that  the  static  pressure  in  the  region  of 
a  nacelle  nozzle  exit  is  given  by  the  empirical  equation 


Cp  =  -0.01  {(uj/u0)  -l)  ,  (20) 

where  u  is  the  free  stream  velocity  and  u  is  the  jet  velocity. 

J 

This  result  is  based  upon  experimental  observations  made  with  a  finite  free- 
stream  velocity  u0  ,  and  cannot,  of  course,  be  applied  to  the  static  case  of 
uo  =  0. 

In  Reference  51,  Payne  was  able  to  derive  this  result  theoretically  for  the  two- 
dimensional  case,  using  the  familiar  analysis  for  potential  flow  constrained  by 
an  obtuse  angle,  as  shown  in  Figure  4. 
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The  analogy  with  an  in  viscid  jet  issuing  from  a  body  is  obvious  from  Figure 
4(a).  When  entrainment  occurs,  some  of  the  free-stream  air  enters  the  jet, 
so  that  the  angle  through  which  it  is  turned  at  the  trailing  edge  is  a  little  less, 
by  the  amount  6  as  noted  below. 


Figure  4.  A  Potential  Flow  Model  of  Jet  Entrainment. 


It  is  easy  to  show  that 


6  =[dn/2(x/t)  )  Uj, 


(21) 


where  n  =  edtrained  air  mass  flow 

jet  mass  flow 


From  this.it  follows  that  the  pressure  decrement  over  the  rear  of  the  body  is 


Cp  -  -K{  (u-/u0)  -l)  , 


where  K 


,08/v[l  +  t'e/n;)2  if  3n/3 (x/t)  =  .08  . 


For 

0  -•  0° 

20° 

40° 

K  =  .0854 

.0206 

.0170 

6(f 

.0143  . 


(22) 
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Kuchemann's  experimental  observation  of  K  =  .01  was  for  axisymmetric 
flow,  whereas  the  above  values  are  for  two-dimensional  flow,  which  is  bound 
to  give  somewhat  higher  suction.  Thus,  the  theoretical  result  can  be  consid¬ 
ered  reasonable,  without  including  boundary  layer  effects.  Also,  since  it  is 
based  on  the  Rouse  et  al  (Reference  37)  measurement  of  3 n/1(x/t)=.08  t  which 
strictly  applies  only  at  9  =  90°  ,  the  theoretical  result  probably  has  errors 
attributable  to  the  change  of  the  mixing  derivative  with  wall  angle. 

A  more  dramatic  example  of  the  potentially  large  effect  due  to  mixing  occurs 
when  a  jet  issues  from  a  bluff  body,  as  illustrated  in  Figure  5. 


Figure  5.  Jet  Issuing  From  the  Stern  of  a 
Bluff  Body. 


In  Reference  30,  Payne  and  Anthony  show,  both  theoretically  and  experimentally, 
that  viscous  mixing  effects  in  this  case  are  so  large  that  the  drag  coefficient  of 
the  body  is  increased  by  several  hundred  percent  when  the  jet  is  operating.  Viewed 
alternatively,  half  the  jet  thrust  is  lost. 
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/mother  effect  of  mixing,  and  in  some  ways  the  most  important,  is  the  effect 
which  it  has  on  the  jet.  Until  the  advent  of  jet  flaps  and  ground  effect  machines, 
aerodynamics  was  primarily  concerned  with  the  reaction  force  obtained  from  a 
jet,  or  its  momentum  flux,  neither  of  which  was  significantly  influenced  by  mix¬ 
ing.  In  fact,  one  of  the  fundamental  theorems  of  mixing  is  that  momentum  is 
conserved  if  the  process  takes  place  at  constant  pressure.  Total  pressure  is 
not  conserved,  however,  so  that  when  we  are  concerned  with  jet  characteristics 
some  way  downstream  of  its  nozzle,  the  conditions  bear  little  relationship  to  the 
same  problem  in  inviscid  flow. 

This  observation  leads  to  a  further  aspect  of  viscous  mixing :  the  apparent  ’’skin 
friction"  loss  involved  in  fluid  shear. 

Skin  friction  results  in  a  drag  force  which  is  defined  as 


DSF  ~  cf  *  *pUo  SWET  , 


(23) 


*«* 


\ 


a  typical  value  for  Cf  being  .005,  for  a  smooth  solid  surface  in  air,  although 
it  varies  with  Reynolds  number,  of  course. 

We  can  obtain  an  equivalent  value  for  Cf  when  air  flows  into  contact  with  an 
air  surface  (Figure  6)  rather  than  a  solid  surface.  At  constant  static  pressure, 
the  total  momentum  flux  will  be  constant.  That  is, 


Figure  6.  Velocity  Shear  at  the  Edge  of  a  Jet. 


fc. 
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The  faster  stream  of  air  has  given  up  energy  in  order  to  achieve  this  result, 
however. 

For  convenience,  and  because  it  greatly  simplifies  the  equations,  a  one-dimen¬ 
sional  analysis  is  used  to  determine  the  value  of  Cy  .  If  the  faster  stream  of 
air  is  regarded  as  a  two-dimensional  jet  of  thickness  t  ,  and  if  we  define  a 
ratio 


entrained  air  mass  flow  (25) 

n  ~  - - - - - , 

initial  jet  mass  flow 

then  the  air  mass  in  the  mixing  zone  is 

(n  +  1)  mQ  . 

For  conservation  of  momentum 

'mouo+™oul  *  ™0  (n+1)  “m  • 

um  =(uQ  +  nuj  )/(n  -hi)  .  (26) 


The  loss  of  momentum  in  the  original  jet  air  mass  flow  is 


DSF  =  Vo  -  'V*  : 

therefore, 

Cf  =  2  (t/x)  [n/(n+l)){l-(u2/uo)]  .  (27) 

For  a  two-dimensional  jet,  in  the  zone  of  flow  establishment  (Reference  37  for 
example)  , 

n  «  3n/3 (x/t)  .  x/t  ; 

therefore, 

Cf  =  2(3n/3(z/t))  .(l-(u2/u0))/(l+(x/t)On/3{x/t}))  .  (28) 


«N» 

t  ■■■■>. 
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Reuse  et  al  (Reference  37)  measurements,  modified  (Reference  51)  for  velocity, 
give 

3n/a (x/t)  =  .08  (l  -  (Uj/Uq))  ; 

therefore, 


Cj>  -  0.16  based  on  relative  velocity  , 

or  approximately  thirty  times  as  great  as  the  value  of  .  005  obtained  for  air  flow¬ 
ing  over  a  smooth,  solid  boundary. 

This  result  is  significant  and  one  which  is  often  not  sufficiently  appreciated  by 
engineers,  perhaps. 

TRANSFER  OF  POWER  BY  VISCOUS  MIXING 


SECONDARY 

FLOW 


PRIMARY 

FLOW 

0 


Figure  7.  Basic  Geometry  for  Calculating  Power 
Transfer  Between  Two  Streams. 


To  study  power  transfer  from  a  primary  to  a  secondary  flow,  the  imaginary 
line  O  -  O  is  drawn  along  a  jet  streamline,  as  shown  in  Figure  7.  Thus,  no 
flow  crosses  this  line.  From  Equation  (26),  for  constant  static  pressure  Apj , 


urr/uj  =[l+n  (uj/i^)Y (n+1)  =[l+n2  (Aj/Aj)]/(n+l))  -  (29) 
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If  A/?  is  the  total  pressure  rise  which  occurs  in  the  secondary  flow,  the  power 
added  to  it  is 


nhj  Atf/p  .  (30) 

For  constant  static  pressure, 

=(%p u2-/(n*D2)[l  -  (uj/u  j)  ( 1  +  (u1/u.)  {1  *  2n })  .  (31) 

4  V  *  J  1  J 

Thus,  the  secondary  flow  power  increase  is 

tPj  * (!fn/hMj/(n+l)2)(l-(u2/Uj)}  (l+fuj/uj)  { I+2n})  .  (32) 

The  total  pressure  drop  in  the  primary  is 

AHj  -  W.  iJ-(urr/uj)2) 

:s[3fpuy(n-fl)2)n[l-(u2/uj))[l-t-(uj/u;j)  {l+2n})  .  (33) 

The  power  lost  by  the  primary  fluid  is  therefore 

APj  -[kmjUytn+D^n  [l-(uj/uj))  [n{l+(u2/uj)}+2)  .  (34) 

Thus,  the  efficiency  of  power  transfer  is  obtained  by  dividing  Equation  (34)  into 
Equation  (32).  That  is, 

n  =[l+(u2/Uj){l+2n})/[2+n  {l+fUj/vj}})  .  (35) 
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The  entrainment  ratio  n  is  defined  by  the  mixing  length  .  For  the  two- 
dimensional  case,  from  Reference  37,  in  the  zone  of  flow  establishment, 


(36) 


Equations  (35)  and  (36)  permit  the  efficiency  of  power  transference  to  be  plotted 
as  a  function  of  Uj/uj  and  l^/tj  .  The  secondary  air-mass  flow  is  nhj  , 
of  course,  and  the  total  secondary  flow  power  is  given  by  Equation  (32).  More 
conveniently,  writing 

rn-  =  p Am.  =  pbtu  .  , 

J  J  J  j 


Pj/l^boul  =  (.04{l-(u2/Uj)}2 /(n+1)2)  [l+(ui/uj){l+2n))  ,  (37) 


where  b  is  the  width  of  the  mixing  zone  normal  to  the  plane  of  the  paper. 


The  ratio  is  a  fixed  quantity.  From  Reference  37,  an  average  value  is 


Thus,  all  the  important  variables  are  a  function  of  ul/u0  only. 
The  specific  mass  flow  is  obtained  from  the  relationship 


rvn./lhpu.  =  .08  [l-fu./u  .)) 
J  tJ  <7 


(38) 


Efficiency  is  plotted  in  Figure  8  and  the  specific  power  loading  in  Figure  9.  Even 
at  low  values  of  the  velocity  ratio  uj/uj>  the  efficiency  of  the  transfer  process 
is  quite  reasonable.  It  is  evident  that  an  explanation  of  the  very  low  efficiencies 
currently  measured  for  eductors  must  be  sought  elsewhere. 
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THE  CHARACTERISTICS  OF  JET  FLOW 


EXPERIMENTAL  OBSERVATIONS 


Before  discussing  the  effects  of  mixing  on  eductor  performance,  it  would  be  well 
to  review  a  few  of  the  known  facts  about  the  phenomenon  of  jet  flow. 

As  an  initial  generalization,  jet  flow  may  be  divided  roughly  into  four  separate 
categories,  the  Reynolds  number  (Re)  of  a  jet  determining  the  category  which 
best  describes  its  behavior. 

Reynolds  number  is  defined  as 


Re  -  ut/v 


(39) 


where  u  =  the  jet  velocity. 

v  =  the  kinematic  viscosity  u/p  . 

U  =  the  viscosity  of  the  fluid. 

p  =  the  density  of  the  fluid. 

t  =  a  characteristic  length,  usually  the 
minimum  dimension  across  the  jet 
nozzle. 

At  very  low  velocities,  the  flow  of  a  jet  is  laminar  and  is  described  by  the 
Navier- Stokes  equations,  which,  for  two-dimensional  flow,  are 


3u/3£  +  u(du/'dx)+  v(%u/ty)  =  -2/pCdp/Zx)  +(\l/ pjv2  u 
dv/dt  +  u('bv/'hx)+  v(<)v/9 y)  =  -l/pCbp/lx)  +( \i/ p)V2  V 


(40) 


When  the  Reynolds  number  is  very  small  (  Re  <  1  )ithe  inertial  terms  on  the 
left  hand  side  can  be  neglected,  and  the  simplified  equations  describe  ’’creeping 


flow. "  In  this  region, theory  gives  good  agreement  with  experimental  observations. 
Needless  to  say,  this  is  not  a  flow  condition  which  has  much  practical  application. 

Below  Reynolds  numbers  in  the  range  25  -  1000,  the  flow  pattern  of  a  jet  is  char¬ 
acteristically  laminar,  and  of  the  type  illustrated  in  Figure  10. 


(b)  Laminar  Jet  From  a  Thin  Tube. 


Figure  10.  Typical  Laminar  Flow  Jets. 


A  few  solutions  exist  for  idealized  laminar  jets,  but  again  the  Reynolds  number 
is  so  low  that  the  results  could  be  applied  only  to  small  models,  since  He  -  25  to 
1000  implies  0.004  <ui  <  0. 16  ft2/sec.  at  sea  level.  The  usefulness  of  laminar 
solutions  is  extended  by  three  considerations,  however.  At  high  Reynolds  numbers, 
when  the  jet  is  turbulent,  it  is  surrounded  by  a  laminar  "sheath"  in  the  region  where 
the  local  velocity  is  less  than  half  the  maximum  jet  velocity.  This  is  presumably 
because  of  the  low  velocity  gradients  involved,  away  from  the  actual  shear  region. 
Secondly,  it  is  theoretically  possible  to  assume  an  effective  "eddy-viscosity"  co¬ 
efficient  which  permits  turbulent  flow  to  be  represented  by  the  Navier-Stokes  equa¬ 
tions.  The  "eddy-viscosity"  has  a  much  larger  value  than  the  true  (laminar  flow) 
viscosity,  and  is  not  a  linear  function  of  shear  velocity  gradient,  so  that  its  use¬ 
fulness  is  limited  at  the  present  stage  of  theoretical  development.  Thii’dly,  it  has 
been  suggested  that  although  the  near  field  mixing  process  is  basically  turbulent, 
the  jet  may  be  laminar  for  a  short  distance  downstream  from  its  nozzle.  The  laminar 
theory  may  be  applied  in  this  area. 

Returning  to  the  description  of  jet  characteristics  as  a  function  of  Reynolds  num¬ 
ber,  it  has  been  established  that  the  creeping  flow  which  occurs  at  very  low  values 
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of  Re  becomes  laminar  as  Re  increases.  A  further  increase  to  the  region 
100  <  Be  <  6000  causes  the  jet  to  become  "ragged"  and  periodic,  as  illustrated 
in  Figure  XI. 


Figure  11.  Periodic  Jet  Structure  in  the 

"Transition"  Reynolds  number 
Region. 


Since  the  transition  Reynolds  number  region  corresponds  to  the  range  0, 016 
<  ut  <  o.  375,  small  models  could  experience  this  phenomenon.  If  the  jet  noz¬ 
zle  of  a  peripheral  jet  GEM  model  were  one-eighth  of  an  inch  wide,  for  example, 
the  jet  velocity  would  have  to  be  about  36. 0  ft/sec.  to  insure  that  periodic  jet 
flow  did  not  occur.  This  corresponds  to  a  cushion  pressure  of  1. 0  -  1. 5  Ib/ft.  . 
Thus,  lightly  loaded  GEM  models  could  experience  scale  effects  which  would 
make  their  performance  quite  different  from  that  of  geometrically  similar  but 
larger  models. 

The  "super-critical"  Reynolds  number  region  starts  at  roughly  Re  -  6000, 
and  is  the  region  of  interest  in  most  practical  applications.  For  a  two-dimensional 
jet,  the  flow  picture  is  indicated  in  Figure  12.  It  should  be  noted  that  there  are  two 
distinct  flow  regimes  and  that  the  theoretical  analysis  is  quite  different  in  each 
regime. 

THEORETICAL  ANALYSIS  OF  TURBULENT  JET  FLOW 


At  present,  the  theory  of  turbulent  jet  flow,  or  indeed  of  any  turbulent  mixing 
process,  is  rather  tenuous  and  unsatisfactory  because  of  the  mathematical  diffi¬ 
culties  introduced  by  the  turbulence.  However,  it  may  be  useful  to  briefly  re- 
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view  some  approaches  which  have  been  taken  to  see  how  experimental  data  can 
be  correlated  and  analyzed. 

For  flows  which  are  nearly  parallel  to  the  x  axis  (  y  an  order  of  magnitude 
less  than  u  ),  the  Navier-Stokes  equations  can  be  simplified  by  assuming  quasi¬ 
static  flow  ( 3u/  %t=0  )  and  by  neglecting  the  normal  velocity  component  (  y  ) 
entirely.  This  gives  the  approximate  boundary  layer  equations  of  Prandtl, 


u  (3u/3 x)  +  v  Cbu/ty)  =  -(1/p)  (dp/dx)  +  (p/p)  (3 2u/2x2) 


3u/3x  +  3y/3 y  =  0  . 


Figure  12.  Approximate  Flow  Picture  for  a  Turbulent  Jet. 


N 


The  static  pressure  p  is  a  function  of  the  external  flow  field,  and  is  therefore 
assumed  to  bo  known.  In  the  most  usual  case,  dp/dx  =  o,  and  it  then  follows 
from  Equation  (41)  that  momentum  is  conserved.  That  is, 

■h* 


j  u2ay  =  const. 


(42) 
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Equation  (41)  has  been  solved  for  laminar  fiow  by  many  investigators,  for  both 
free  and  fixed  boundary  problems.  It  can  only  be  applied  to  turbulent  problems 
by  using  the  fictitious  scalar  ‘'eddy-viscosity'*  coefficient  (  e  )  in  place  of  the 
Newtonian  coefficient  (  v/o  )  for  laminar  flow.  Although  (  v/o  )  is  constant 
in  Newtonian  fluids,  (  f-  )  is  not,  so  that  the  resulting  analysis  is  consider¬ 
ably  in  error.  Also,  the  actual  value  of  the  eddy- viscosity  can  only  be  defined 
as  the  value  which  gives  the  right  answer!  Thus,  it  performs  the  function  of  an 
empirical  constant. 

Other  and  quite  different  approaches  have  been  used,  employing  either  dimen¬ 
sional  analysis  or  momentum  theory  analysis.  In.  each  approach,  however, 
empirical  constants  appear  and  must  be  evaluated  by  appeal  to  experiment. 

Perhaps  the  most  elegant  analysis  is  based  upon  momentum  considerations.  An 
excellent  example  is  that  presented  by  Holdhusen  in  the  discussion  section  of  Ref¬ 
erence  37.  Holdhusen  points  out  that  there  are  three  basic  assumptions  needed 
to  derive  the  momentum  equations:  hydrostatically  distributed  pressure,  con¬ 
stancy  of  momentum, and  the  dynamic  similarity  of  the  diffusion  process  at  all 
positions  along  the  jet.  The  first  assumption  is  easy  to  justify,  the  pressure  dif¬ 
ferences  arising  as  a  result  of  flow  in  the  y  direction  being  less  than  1  per¬ 
cent  of  the  jet  dynamic  pressure.  The  second  assumption,  constant  momentum 
flux ,  follows  from  the  assumption  of  constant  static  pressure  and  Newton's  sec¬ 
ond  law.  Alternatively,  It  may  be  derived  from  the  Prandtl  equations,  as  men¬ 
tioned  earlier.  The  third  assumption  is  simply  an  assumption,  and  although 
widely  used  in  this  and  other  problems  (such  as  the  wall  jet  problem,  for  exam¬ 
ple)  it  can  only  bo  justified  by  appeal  to  experiment. 

Defining  n  =  y/x  ,  (43) 

the  "similarity"  hypothesis  implies 


u/umax  =  f  • 

In  the  zone  of  established 


P  tuo 


p  U' 


max 


(44) 

flow,  conservation  of  momentum  gives  the  equality 

f2  (w)  dy  ; 


therefore, 


■b* 


(r/t)j  f2( r\)  dy  «  ^/‘W’2  . 


(45) 
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Because  of  the  similarity  hypothesis,  therefore. 


SioaAo  =  (t/x)if  /  (/  f2(r])  A*  •  (46) 

For  the  zone  of  flow  establishment,  in  order  to  use  the  similarity  hypothesis,  we 
must  take  the  origin  y  =  0  to  be  on  the  edge  of  the  diffusion  zone.  Then  for 
conservation  of  momentum  in  this  region,  where  umax  ~  uo  , 

oo  Q 

P tu1  =  2 j  pu2  f2  (r))  dy  +  2 j  p u2  dy  ; 

0  o  0  -y  e.l. 

therefore, 

j-co 

(x/t)\  f2  (n)  dn  +  2y  .  /t  =  1.0  •  (47) 


Thus,  the  width  of  the  potential  zone  is 

+OD 

2y  /t  =  1  -  (x/t)  f  fz  (n)  dr)  ,  (48) 

o. t. 


showing  it  to  be  bounded  by  planes.  The  value  of  x  at  which  it  vanishes  is 
given  by 


+00 

x0/t  =  1/J  f2  (ry)  dr,  .  (49) 

—  OO 

Thus,  knowing  X<?A  ,  we  know  the  value  of  the  momentum  integral.  For 
the  average  value  of 

x/t  =  5.2 

given  by  Rouse  et  al  (Reference  37) , 

+oo 

/  f2  (n)  dr)  =1/5.2 


Baines,  in  the  discussion  following  Reference  37, reported  values  of  xc>/  ^  in 
the  range  5.0  -  7.0  for  circular  jets,  the  value  varying  with  Reynolds  number. 
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Thus,  the  function  /  (  n  )  would  seem  to  be  a  function  of  Reynolds  number,  a 
result  which  is  not  too  surprising  perhaps,  but  which  gives  an  indication  of  the 
difficulties  inherent  in  analyzing  turbulent  processes. 

By  substituting  Equation  (49)  into  Equation  (46), we  have,  for  the  zone  of  estab¬ 
lished  flow , 


umax^uo  (xq/Xj ^ 


(50) 


The  jet  mass  flow  is  obviously 

+°° 


■f® 


m 


<r 


]/ho  =(/  P“  ifMdyVptUo)  =  (uma/u0)  (x/t)  /  f  (r])dr] 


=  (x/t)  I  f  (x\)dr\ 


(51) 


in  the  zone  of  established  flow. 

In  the  zone  of  establishment, 

•f®  -f® 

PU0f( rt)dy  +  2y  q  ^  pii0)/p tuQ  =  x/tf  f(r\)dr ) 

_O0  —00 

+  (2ya.i.  /t} . 

or  substituting  Equation  (48)  for  $/  c  j  , 

-f-CO  «/-00 


=  1  +  (x/t)  (/  f  (r,)  dr\  -  j  f2(r\)  c?n)  . 

—  OO  —00 


(52) 


(53) 


Momentum  analysis  reveals  a  great  deal  about  the  entrainment  process.  How¬ 
ever,  experiment  is  still  needed  to  determine  "arbitrary  constants.  "  In  the  ap¬ 
proach  described  earlier  using  the  Prandtl  equations  .only  one  "constant",  the 
"eddy-viscosity"  (  c  ),was  apparently  needed.  However,  the  assumption  that 
it  was  constant  also  implied  assumption  of  a  velocity  distributibn  which  would  satis¬ 
fy  Prandtl’s  equations;  therefore,  really  two  assumptions  were  made. 


The  same  remarks  are  true  of  all  other  analyses  of  this  type,  such  as  Tollmein's 
(Reference  35)  classical  analysis,  for  examr.le,  which  uses  the  Prandtl  "mixing 
length"  hypothesis  for  defining  viscosity. 


The  case  of  momentum  analysis  is  no  worse, since  the  momentum  integral 
(  if2  (n)  can  be  determined  from  measurements  of  (  xQ  ),and  the  re¬ 
maining  integrals  follow  when  the  velocity  distribution  f  (  n  )  is  defined. 
Many  investigators  have  used  the  normal  error  function, 


u/u 


max 


-y 2/Sn^ 
e 


(54) 


for  analysis  and  for  correlation  of  test  data.  As  pointed  out  by  Reichardt  (Reference 
40),  however,  it  has  the  inherent  weakness  of  assuming  that  the  turbulence  gen¬ 
erated  by  mixing  extends  to  infinity. 

As  a  matter  of  fact,  there  is  no  need  to  make  any  assumption  in  regard  to 
velocity  distribution,  since  it  can  be  measured  quite  accurately  at  various  sta¬ 
tions,  and  non-dimensionallv  plotted  as  a  single  curve,  permitting  the  integrals 
of  /  (  n  )  to  be  evaluated  graphically.  It  is  even  possible  that  the  values  of 
these  integrals  of  f  (  n  )  could  be  used  to  deduce  characteristics  of  the  function 
itself,  and  hence  reveal  something  of  the  mixing  process.  This  is  unnecessary 
from  an  engineering  point  of  view,  however,  since  only  predictive  ability  is  re¬ 
quired. 

The  measurements  of  Rouse  et  al  (Reference  37)  give  the  following  equations 
for  the  two-dimensional  mixing  process  when  suitably  analyzed  to  give  the  re¬ 
quired  "empirical  constant, "  and  using  the  Gaussian  velocity  distribution. 


m-/m  -  1  +  0.D8  x/t 

v  ts 

(x  <  XQ) 

(55) 

=  0.62  (x/t) ' 

(x  >  XQ) 

(56) 

If  the  mas 3  of  entrained  air  is  defined  as 


so  that 
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m  ./m  =  (n  +  1)  or  n  =  [(m./mn)  -  2)  , 
jo  j  u 

then  n  =  .08  x/t  (x  <  xQ)  (57) 

=  f 0.82  (x/t)*  -  j)  (x  >  x  )  , 
v  o 

and  the  entrainment  per  unit  (  x/t  ) 

dn/d(x/t)  =  .08  (x  <  xQ) 

k 

=  0.31/ (x/t r  (x  >  xQ)  .  (58) 


These  results  are  plotted  in  Figures  13  and  14.  It  is  interesting  to  note  that  the 
value  of  3«/3  (x/t)  is  constant  for  x/t  <  5.2,  which  is  the  region  at  great¬ 
est  interest  for  many  practical  applications,  although  not  for  the  eductor  problem. 


Figure  13.  Variation  of  Entrainment  Ratio 
n  with  x/t  . 


33 


I 


Figure  14.  Variation  of  Entrainment  Function 
Derivative  3n/3  (x/t)  with  x/t  . 


When  a  jet  exhausts  parallel  to  an  axial  free  stream  ( uQ  ), following  a  suggestion 
made  by  Kuchemann  in  Reference  83,  it  can  be  assumed  that  the  rate  of  mixing 
is  proportional  to  the  shear  velocity  (  u~uo ).  This  amounts  to  multiplying  the 
x-ordinates  by  the  factor 

uy(uj-u0)  *  l/[l  -  (%/Uj)}  (59) 


so  that  the  entrainment  functions  become,  from  Equations  (57)  and  (58), 

n  =  .08  (x/t)  (l  -  (u0/Uj))  (x  <  x0) 

=  0.S2  {(x/t)  (1  -  {Uq/u.}))3*  -  1  (x  >  xQ) 


(60) 
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(61) 


dn/d(x/t)  =  .08  (2  -  ('n/Wj*')  (s  <  xQ) 
=  {0.31/  (x/t)H\l  -  (u/u.)) 


(x  >  XQ)  . 


When  a  jet  exhausts  normal  to  a  free-stream  flow,  a  transverse  flow  field,  as 
such,  cannot  occur  across  it.  Instead,  an  initially  normal  airflow  (of  velocity 
u0 )  will  be  arrested  by  the  jet, resulting  in  a  local  static  pressure  differential 
of  AP  -  ho ,  or  will  be  deflected  by  the  jet  in  a  manner  which  will  depend 
upon  the  strength  and  size  of  the  jet  and  the  boundary  conditions. 
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Chapter  5 

THE  EFFECT  OF  MIXING  STATIC  PRESSURE 


A  basic  phenomenon  is  the  effect  of  mixing  which  occurs  at  a  static  pressure 
which  is  not  equal  to  ambient. 


Figure  15.  Localized  Mixing  in  a  Jet. 


At  constant  static  pressure  (  AP]),two  streams  of  fluid  will  mix  in  accordance 
with  the  law  of  conservation  of  momentum,  as  discussed  previously.  Assume 
now  that  a  jet  entrains  a  secondary  flow  within  a  control  volume  (  1  -  m),  as 
shown  in  Figure  15. 

If  the  static  pressure  in  the  control  volume  is  tyj  ,  if  bj  is  the  jet  mass 
flow  in st  station  (l),arri  if  {  2+n )  is  the  mass  flow  out  at  (m),  then  from  Equa¬ 
tion  (29) , 

-(1  +  n  ( uj/uji)/(n  +  1)  .  (62) 


So  long  as  the  mixing  pressure  continues  at  Apj ,  there  will  be  no  change  of 
momentum.  Eventually  the  jet  must  return  to  ambient  static  pressure  (  Ap  =  o), 
however,  and  if  this  is  achieved  without  loss  of  total  head, 
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A Pj  +  ipu£  =■  Hpu\ 


Thus,  the  final  (  u2  )  velocity  is 

u2  =  (2/p  (bp1  +  bou2))*  t  (63) 

and  the  final  momentum  flux 

J2  =  (n  +  1)  ^ju2 

•  (n  +  1)  m.  ((2/p)  5p2  +  .  (64) 

Expressing  this  as  a  ratio  of  the  initial  jet  momentum  flux, 

J2/Jo  m  +  V  (“P,  +  (um/u4)2)i  ■  <65> 

i  J 

where  Lp^  =  bp^/hpii2  . 

Substituting  Equation  (62)  for  u^/u^  , 

J2/J0  =  ((n  +  l)2  t$1  +  (1  +  n  {uj/uj})2}*  .  (66) 

This  equation  offers  some  interesting  special  cases,,  For  an  ideal  static  aug- 
mentor  with  no  internal  losses, 

HPU2  -  -  A p1  ; 
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-  4. 


therefore, 

if 

Uj/Uj  =  (-bp2) 

L  L 

and  J2/Jo  m  ((n  +  2)2  *  {I  +  n  (~LPj} *2)  *  (67> 


This  function  is  plotted  in  Figure  16  for  n  =  10. 0*  It  is  evident  that  there  is 
an  optimum  value  of  Ap^  for  maximum  increase  in  momentum  flux.  Static 
eductors  should  obviously  be  designed  to  take  advantage  of  this  optimum. 


For  static  entrainment  in  a  jet . 


u  ^  —  0 


therefore, 

J  /J  —  (n  +  1)  +  1 

2  o  i 


That  is,  J 2^0  >  1.0  for  greater  thru  ambient  mixing  pressure,  the  reverse 
of  a  static  augmentor. 

An  example  of  such  entrainment  occurs  on  the  high-pressure  side  of  a  jet  flap, 
which  is  deflected  through  a  large  angle.  Close  to  the  nozzle  the  relative  static 
pressure  may  be  taken  as  the  free-«tream  dynamic  head  %p u2  on  the  high- 
pressure  side  of  the  jet.  On  the  rear  side  of  the  jet,  the  static  pressure  will 
be  ambi-nt,  or  somewhat  lower.  Thus,  the  mean  effective  jet  static  pressure 
will  be  between  these  two  extremes,  and  will, in  fact, be  the  arithmetric  mean  if 
the  jet  is  thin  in  relation  to  its  local  radius  of  curvature.  Thus,  approximately 

Lp2  «  Ap  .  =  idfpu2)  ; 

therefore, 

api  =  *<w2 


and  J2/J0  =  (*  (n  +  l)2  Cu^'uj)2  +  l)h  .  (69) 

This  equation  does  not  account  for  the  momentum  lost  by  the  air  which  is  entrained 
by  the  jet,  a  drag  force  which  is  .  When  this  is  included , 


AUGMENTATION 


=  [h(n  +  l)2(u°/uj)2  +  l)*  -  n  (u^/uj)  . 


(70) 


yj 


NET 


Note  that  as  u^u.  +  O  f 


JJJ0\  +  1  -  n(u  /u  .) 

2  °\net  0  0 


As  n  ■*  0 


(71) 


J2^Jo 


NET 


(%  (uc/uj)2+  2)‘  , 


(72) 


which  describes  the  effect  of  the  pressure  energy  in  the  jet. 

In  general,  the  ratio  J2^o  will  be  less  than  unity.  Thus,  this  mixing  loss 
is  the  reason  for  the  so-called  '•thrust  loss  anomaly"  of  the  jet  flap. 

While  it  is  not  the  primary  purpose  of  this  report  to  discuss  jet  dra*  eful 

to  note  that  the  variation  of  (  n  )  with  uq/uq  can  obviously  be  accoiu*. _ .<*i  using 

results  of  the  type  given  in  the  foregoing  chapter  on  jet  characteristics. 

Free-air  entrainment  is  by  no  means  the  only  cause  of  "jet  drag. "  Apart  from 
"apparent"  force  changes  due  to  the  change  of  local  static  pressure  near  the  jet 
nozzle,  diffusion  losses  occur  on  the  pressure  side  of  a  nozzle  which  has  a  pres¬ 
sure  differential  across  it.  Also,  in  some  configurations,  the  jet  can  increase 
the  size  of  the  wake.  During  recent  experiments,  for  example,  a  jet-flapped  wing 
generated  a  wake  three  times  as  deep  as  the  wing  wake  with  the  jet  inoperative. 

Another  interesting  form  for  Equation  (66)  is  appropriate  to  a  jet  issuing  parallel 
to  a  free-stream  flow.  If  the  free-stream  velocity  is  uQ  ,  then  for  zero  losses 

*P  U2  =  &p2  +  %Quj  ; 
therefore, 

(ui/uj)2  -  <VV2 " 

and  ~  (^n  +  +  ^  +  n  ^U</U^Z  ~  ^Pl***)2)**  -  nfug/uj).  (73) 
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The  thrust  ratio  is  always  less  than  unity  if  the  mixing  pressure  is  greater  than 
ambient,  but  augmentation  can  occur  if  Ap,  <  0.  As  indicated  in  Figure  17, 
the  maximum  attainable  augmentation  diminishes  with  increasing  forward  speed. 

When  &p2  =  0,  =  1.0  always. 

Since  the  static  pressure  behind  a  wing  or  nacelle  is  usually  greater  than  am¬ 
bient,  Equation  (73)  describes  a  "jet  drag"  which  is  the  same  as  that  discussed 
in  the  beginning  of  this  report.  Instead  of  using  a  potential  flow  analysis,  a  dif¬ 
ferent  approach  has  been  used  to  obtain  the  same  result. 


Figure  17.  Effect  of  Forward  Speed  on  Augmentation 
for  n  =  10. 
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t> 

Chapter  6 

GENERAL  ONE-DIMENSIONAL  FLOW  THEORY  OF  AN  EDUCTOR 


TEE  BASIC  EQUATIONS 

As  already  discussed,  the  ability  of  an  eductor  to  increase  the  primary  jet  thrust 
is  predicated  upon  the  mining  pressur£'s(pj  )  being  less  than  ambient,  so  that  dif¬ 
fusion  must  take  place  in  order  to  bring  the  static  pressure  of  the  fully  mixed  jet 
up  to  the  exit  pressure.  There  is  an  optimum  value  for  the  mixing  pressure  ( pj  ) 
which  gives  maximum  augmentation  ,  and  this  leads  to  the  concept  of  an  "optimum 
augmentor. "  This  is  the  augmente  d  which,  for  a  given  diffuser  efficiency  , 
gives  the  absolute  maximum  augmentation  possible.  The  ability  to  generate 
equations  for  this  optimum  reduces  the  computational  load  by  an  order  cf  mag¬ 
nitude  for  a  static  eductor.  When  the  upstream  and  diffuser  exit  pressures  are 
not  ambient,  the  saving  is  even  greater, of  course. 


Let  J 


momentum  flux 


✓* 

m 


a  mass  flow 


P  -  a  total  pressure 

p  =  a  static  pressure 

n  =  a  ratio  of  entrained  to  primary  mass  flow 


The  suffixes  j  ,  1,  2,  and  m  refer  to  the  respective  positions  in  Figure  1. 
The  suffix  a  refers  to  ambient  conditions,  as  when  the  primary  exhausts 
to  ambient. 


The  diffuser  efficiency  is  defined  as 


r\D  =(<p2  V-  *p up  -  Pj)/%pw£ 


(74) 


so  that  it  experiences  a  dynamic  loss  (l-nD)  ipu*  ,  and 


hm\  =  P2  -  P2  *  "d  • 


(75) 
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Substituting  Equation  (62)  for  u 


m  * 

Hqu\  =  pj  -  p2  +  r\D  %pu2  (l  +  n  (Uj/Uj)}2/  (n  +  l)2  •  (76) 

To  non-dimensionalize ,  for  convenience,  divide  by  %(>u2.  ,  so  that 

v 

P„  m 

and  Equation  (76)  becomes 

(u2/'uj)2  ^  Pi  ~  $2  +  nD  i1  +  n  (“l/uj))2/  (n  +  1)2  •  (77) 

From  Figure  1 , 

po  =  Pi  +  (ui/uj)2  ; 

therefore, 

u2/uj  *  (P  -  p2J*  . 

Substituting  into  Equation  (77)  , 

Ug/uj  =  [(pj-Pp)  +{T\D/(n+l)2}{l+n  (Pq-P^)2)^  .  (78) 

The  augmentation  ratio  is 

jp/Ja  «  (n+1)fhju2^jaUja  **  (™/™ja)2(u2/uj)  'n  +  1)m  (79) 

Substituting  Elation  (78)  for  Wg/uj  , 

=  (fn/fnja)2  i(Pi~P2)  (n  +  1)2  +  nD  {1  +  n  (~P c~p )2f  ,  (80) 
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Since  (njMja)2 


buyip.-pj 


(uju .  )2 
i  3d 


t 


(81) 


if 


AP. 

3 


» j( (pj~P2^  ( n  *  ^  *  n 


fpe-p3;*)2)%  . 


(82) 


Equation  (82)  is  a  complete  (one-dimensional  flow)  statement  of  augmentor  per¬ 
formance,  although  the  independent  variables  Pj  and  n  are  rather  incon¬ 
venient  in  this  form. 

OPTIMUM  MIXING  PRESSURE 


The  optimum  mixing  pressure  parameter  is  that  value  of  Pj  which  gives  the 
maximum  value  of  J ^/^with  PQ  and  p*>  being  fixed.  Differentiating  Equa¬ 
tion  (82)  with  respect  to  pj  and  equating  to  zero, the  optimum  value  of  Pj 
is  found  to  be 


Pi 


1  opt 

Substituting  for  Pj 


*  -  n2ny[!n  +  l)2  -  n\)*  . 

into  Equation  (82)  , 


(83) 


opt 


opt 


[(n  +  1)/LP.){dd  +  (Pa‘P2)i(n  +  l)2-  n1^} y/((n+l)  -n2nD)' 


L  ’  v*. 


opt 


2  +  (fn+V/&P.)2  P  -  ((n+V/&P.)2p  . 


(84) 


opt 

fP? 


r  c2 
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This  is  an  interesting  result.  The  effects  of  upstream  and  diffuser  exit  pressure 
appear  as  simple  Southwell  coefficients.  It  is  obvious  that  an  increase  of  inlet 
total  head  will  increase  ^2^ a*  while  an  increase  in  back-pressure  will  decrease 
it,  as  should  be  expected.  Also,  for  ambient  inlet  total  head, the  thrust  falls  to 
zero  when 


*P2  *  *  1)2  “  *  <85) 

It  is  important  to  note  from  Equation  (83)  that  the  back  pressure  does  not 
influence  the  optimum  value  of  Pj  .  This  means  that  an  optimum  injector  is 
merely  an  optimum  eductor  with  a  larger  diffuser. 

THE  OPTIMUM  STATIC  EDUCTOR 

The  preceding  section  shows  that  the  optimum  eductor  is  directly  related,  by 
appropriate  Southwell  coefficients,  to  the  optimum  static  eductor,  for  which 

JJJa\  =  {(n  +  1)  AP-  {(n  +  l)2  -  n2nD}*)  .  (86) 

I  opt 

In  the  limit  n  »  , 

*  M1  -  V)%  •  <87> 

opt 

and  as  n  0  , 

VJa  *  •  <88> 

opt 

Equation  (86)  is  plotted  in  Figure  18,  where  the  dominating  effect  of  diffuser  ef¬ 
ficiency  is  clearly  seen.  Use  of  boundary  layer  control  or  similar  techniques  to 
give  very  high  diffuser  efficiencies  makes  possible  augmentations  as  high  as  3  or 
4,  or  possibly  even  more. 

Calculation  of  geometrical  variables  which  give  optimum  augmentation  follows: 
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AUGMENTATION  RATIO  4> 


Since 


=  -«2ip  /{(n  +  l)2  -  n2r\D)2  =  Ap^/fAP^-Ap-j)  , 


A Pj 

=  -«2n2  /( 

opt 

Api/AP- 

ii 

'•§’  i 

opt 

opt 


/(1+  Ap 


opt 


=  -n2r\fi/({(.i  +  l)2  -  n2T\D}2  -  n2 r\2) 


(89) 

(90) 

(91) 


This  relationship  is  plotted  in  Figure  19.  The  value  of  AP^.  in  Equation  (86) 
can  now  be  calculated.since 


AP-  =  AP/r AP.  -  A pj  =  1/  (l-fAp  /APJ)  .  (92) 


If  we  substitute  this  in  Equation  (86)  , 


V'a 


opt 


=  [(n  +  l)2  -  n2n£l)3//2  (f>i  -f  1)  )/ 


({fn  +  l)2  -  n2r\D)2  -  n2 np)^> 


which  is,  of  course,  the  equation  actually  used  to  compute  Figure  18. 
Now, 


(93) 


«  m  =  r-Ap:;%=({  -aPi/ap.  }/u  -  ^/apj})*  .  (94) 

The  area  ratio  of  the  mixing  tube  is  obtained  from  the  equation  for  continuity  of 
mass  flow , 


AjAi  =  ((n  +  l)2  An-/A2)/(l  +  n2(A./AJ)) 


(95) 


The  diffuser  area  ratio  *s  f°und  from  the  relationship  of  Equation  (75), 

substituting  so  that 
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*P,  *  ip  Vou|M2/VJ  ”  *0“2 


"d  'VV2  * 3  -  *1  W2  •  <96> 

The  velocity  ratio  can  be  obtained  from  Equation  (78).  More  conveniently,  from 
Equations  (79)  and  (81),  and  (92), 

u../u  =[(fnjm.)2  (n+1  )/(J  /J  ))*=[{l-(bp  /tP.)}  (n  +  l)/JJJa  )  ; 

J  o  ]  J**  2  “  *  J 

therefore, 

k^/Am  «  (n  +  1)  /(n Di(n  -hi)2  -  n2^})*  .  (97) 

The  optimum  area  ratio  values  are  plotted  in  Figures  20  through  22,  and  the 
optimum  equations  are  summarized  in  Table  I. 

EFFICIENCY  OF  THE  OPTIMUM  STATIC  EDUCTOR 

The  total  eductor  efficiency  may  be  calculated  from  the  relationship 

n  =[(J2/'Ja)2  /  (n  + 

=  (J/Ja)z  /  (*  +  V  [1  -  (tPj  /Wj))*  •  (98) 

i 

For  the  case  of  np  -  1.0,  n  +f8/sjf3/4r5!  0. 77  as  «-♦■<»  .  For  any  value 
of  as  n  ■*  o,  n  np  . 

Equation  (98)  is  plotted  in  Figure  23  for  optimum  static  eductors. 
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INLET  VELOCITY  RATIO 


! 


Figure  20.  Variation  of  Optimum  Inlet  Velocity  Ratio 

With  Entrainment  Ratio.  (Note  that  u^/uj  = 
(n)  Aj/Aj). 
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DIFFUSER  EFFICIENCY 


Table  I 

OPTIMUM  STATIC  EDUCTOR  PARAMETERS  AS  A 
FUNCTION  OF  ENTRAINMENT  RATIO  n 


Parameter 


Equation 


Total  momentum 
Actual  jet 
momentum 


lopt  (tm. V2-n%)* 


Total  momentum 
Jet  momentum 
to  ambient 


vjo 


opt 


Limit 

n  o 


Limit 
n  ■*  • 


[r\j/(l-r\D)  )* 


Jet  total  pressure  -  „  .  „  ,  .  _  . 

parameter  ^  "  1-^,/^ 


Mixing  pressure  Ap^/AP .  =-n2ni/( { (n+l)2-n2 n_}2 

^  opt  - 


Intake  area  ratio 


-»2# 


opt 


=(2/nZ))({n+i;2-n2nZ}>2  I/n£ 


+2*2n2f  =  (l/n){(- Ap/AP^/  2-  fcp/AP^.  ))' 


Mixing  tube 
area  ratio 

V4?  j  ’s(l+n)2(A/A1)/(Un?Aj/A1 )  nD 

Diffuser  area 
ratio 

V*m|  vry* 

0  ^  1) 
1/3  UD  -  1) 

0 


"D 


v(yj-y)* 


i 

5 

i 
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Jt 


TOTAL  EDUCTOR  EFFICIENCY  n 


f 


Figure  23.  Total  Efficiency  of  Optimum 
Static  Eductors. 
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DIFFUSER  EFFICIENCY 


DIFFUSER  EFFICIENCY 


Thus  far.  the  fact  that  a  conventional  single-stage  diffuser  can  accommodate  only  a 
velocity  ratio  uo/u^  >  %  has  been  ignored, as  well  as  the  fact  that  diffuser 
technology  is  still  a  crude  and  inefficient  business.  There  are  several  ways  of 
improving  both  efficiency  and  diffusion  ratio;  the  most  obvious  is  boundary  layer 
control  utilizing  the  mixing  section  suction.  To  adequately  discuss  this  and  other 
approaches  is  beyond  the  scope  of  this  report,  however. 

HEATING  THE  PRIMARY  AIR  TO  GIVE  INCREASED  PERFORMANCE 

The  density  (  p  )  of  a  gas  is  given  by  the  equation 


P  «  p/RT 


(99) 


The  momentum  flux  through  the  primary  nozzle  is 


m.u  .{(2/v)(hP.-hpJ^  ,  (10°) 

j  j  J  J  1 

Thus,  if  the  primary  air  is  heated, a  primary  momentum  flux  increase  results  in 
the  ratio 


vyp'fyv  ■  (T/T0>h  ■ 


(101) 


if  the  mass  flow  is  held  constant  by  using  variable  area  primary  nozzles.  This 
increase  can  give  an  important ' 'boost"  capability  to  certain  types  of  eductors. 

v 

When  the  nozzle  area  is  not  variable,  the  mass  flow  varies  inversely  as  (T/TQrt 
so  that  the  primary  momentum  flux  (and  therefore  the  total  thrust)  is  independent 
of  temperature.  This  has  been  well  confirmed  in  a  number  of  experimental  pro¬ 
grams.  Indeed,  so  extensively  have  the  effects  of  temperature  been  investigated, 
one  wonders  whether  the  experimenters  appreciated  the  significance  of  Equation 
(99)! 
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PHYSICAL  IXX3ATI0N  OF  FORCES  ON  AN  EDUCTOR 


The  additional  momentum  flux  of  an  eductor  is  transmitted  to  its  wall  in  the  form 
of  static  pressure.  For  the  optimum  case,  calculations  can  be  made  for  the  rel¬ 
ative  magnitude  of  the  primary  entry  lip,  diffuser  and  mixing  chamber  contribu¬ 
tions. 

Primary  Jet 

Assuming  constant  nozzle  total  bead  , 


?o  -  ^Pj  +  P/lj-xj- 

but  P uZ  -  2(&Pj-&Pj)  ; 

therefore, 

Fq  «  2AjhPj.  (j-  h  {hpj/hPj})  ,  <102) 


where  2Aj CJ?j  is  the  value  obtained  in  exhausting  to  ambient  static  pressure, 
of  course. 

Mixing  Chamber 

Fm  *  -Mj-VAPj 

-  -'•jPPj  (I-  V*i  )  •  <103> 


Diffuser 


-h  ’  -  feul  *  P2>  *  fc“m  +  PlJ  A„ 

For  Ps  -  Pa 


9 


-h  -  -  »A?\  +  (eum  *  Am  ■ 

For  continuity, 

U  A  -  U-4.  ; 
m  m  H  * 

therefore, 

•  <1M) 

Not 3  that  the  diffuser  thrust  is  always  negative. 

Summation  o'  Forces 

ft  is  logical  to  divide  the  component  forces  by  24 ^.APj  ,  the  primary  thrust 
to  ambient. 

Then  Fj/Ja  -  (2-  HUpj/t Pj>)  , 

VJa  '  -*  , 

-  *f W { -1  >-*rW  ’  (105> 


Table  n  summarizes  design  calculations  for  three  eductors  optimised  around 
conventional  conical  diffiisers,  and.  summarizes  the  location  of  the  forces  which 
contribute  to  the  augmentation,  ft  is  notable  that  the  bellmouth  inlet  carries  all 
the  net  thrust  increase,  and  also  most  of  the  "negative  thrust"  generated  in  the 
dlfftiser.  The  importance  of  careful  bellmouth  design  therefore  is  clear. 
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Table  n 

SUMMARY  OF  THREE  "OPTIMUM  EDUCTOR"  DESIGN  CALCULATIONS 


Entrained  air  ratio  » 

5.000 

10.000 

15. 000 

Augmentation  J  2^ a 

2.000 

2.300 

2.600 

Diffuser  efficiency  nr, 

0.960 

0.945 

0.940 

Diffuser  area  ratio  A2^Am 

1.800 

2,200 

2.470 

Diffuser  diameter  ratio  <2 2/^ 

1.340 

1.490 

1.570 

Diffuser  length  ratio 

6.530 

9. 220 

10.920 

Primary  jet  area  ratio  4^/4 . 

14.  370 

31.400 

51.000 

t./Di  (jet  ring  diameter  =£7 ) 

0.017 

0. 008 

0.005 

1 Jt . 
nr,  j 

V»i 

143. 000 

461. 000 

935. 000 

2.500 

3. 680 

4.500 

Mixing  zone  area  ratio  A2/Al 

0.920 

0.920 

0.920 

Mixing  zone  diameter  ratio  d^/dj 

0.  960 

0. 960 

0.960 

Bell  mouth  diameter  ratio  dJJd1  = 

2.00C 

2.000 

2.000 

Mixing  pressure  Ap^/AP^  -0. 190  -0. 150 

Breakdown  of  Location  of  Thrust  Increase 

-0. 110 

Thrust  increase  ( ( J 2/Ja)-l ) 

1.000 

1.300 

1.600 

Primary  jet  nozzle 

+0.092 

+0.053 

+0.036 

Mixing  chamber  wall 

+0.113 

+0. 129 

+ 0 . 138 

Diffuser  wall 

-0.363 

-0. 548 

-0.719 

Bellmouth  lip 

+1. 158 

+1.  358 

+1.  545 

PERFORMANCE  OF  A  STATIC  EDUCTOR  OF  ARBITRARY  GEOMETRY 


In  the  previous  chapter,  the  "optimum  eductor"  was  analyzed  and  defined  as  an 
eductor  whose  geometry  is  optimized  for  maximum  performance.  In  this  section 
the  equations  are  rearranged  to  permit  the  performance  of  a  specific  augmentor 
to  be  calculated  for  static  conditions. 

There  does  not  appear  to  be  an  explicit  solution  for  an  arbitrary  eductor  geometry, 
and  a  solution  must  be  obtained  by  iteration,  or  fcy  plotting  two  functions  against  the 
entrainment  ratio  n  ,  in  order  to  determine  their  intersection. 
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In  an  "optimum  eductor','  both  the  intake  and  exit  areas  are  defined  by  the  entrain¬ 
ment  ratio  (  n  ).  If  the  intake  area  (  Aj  )  is  already  defined,  there  is  still  an 
optimum  mixing  pressure  and,  consequently,  an  optimum  diffuser  area  ratio. 

The  equations  defining  this  optimum  are  defined  later  in  this  section. 


Performance  Equations 


Stoce  • 


ip3  +  n0  tm$  i'VV2  ”  , 


or  *Pj  »  -<u/u.)2  (VW2  -3> 


(106) 


Now  J2/JQ  -  (n  +  1)  (uJuj)  ,  (107) 

(VA/Am)2'^/(1  +  n)2  *  (108) 

But  (J/J0)2  *  AP2  a  (1  *  n  (-& pj)*)  .  (109) 


Substituting  Equation  (108)  for  bp  ^  , 

tJ2/J0)2  (2  +  l»t(VV!'I)  il-(n2n,y'(l+n)2)})  (110) 

-2  V/V  (nnj/d+n))  (WV2*  :f-nD  =  0  ; 

therefore, 

V-'o  !n0-(V‘12-|2>*)/(3-('>2/(n+2;2)  (111) 

<vfw2»  • 

Now 

_  % 

n  =  j^/A  .  (-bpj) 
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(112) 


-(rVV^V  {nD  M/V2-I1%}/rn  + 


from  Elation  (108).  Therefore, 


J2^Jo  =  [(An/A2*  (A/Aft  ln(rt+1)  ~  (An/A2) 

By  evaluating  Equations  (111)  and  (113)  for  a  range  of  n  values, 
coincidence  may  be  determined,  and  hence  both  J2/J0  and  n  . 
obtained  by  the  relations 

Ap.  *=  n2(A  ./A.)2 
1  j  l 


bP.  =  1  +  AjL 
J  1 

j  /ja  =  (j2/j0)  (1^j}  *  w  7  <J+n2^/V2} 


Optimum  Diffuser  Area  Ratio 

From  Equation  (83)  of  the  previous  chapter  , 


&p2 


opt 


»--n2 r)2  /  {(n+l)2-n2r\D): 


and  =  ( (n+l)2-n2r\D)3/2  (n+1)  ( r\D ft  /  {(n+l)2-n2r\v)2- 


From  Equation  (114), 


n  =(A  /A-)(-Ap1  ft 
opt  1  v  opt 

=  (A2/Aj)mD  /{(n+1)2-  n2nD )  . 


(113) 

the  point  of 

VJa  18 

(114) 

(115) 

n2n*  .  (116) 

(117) 
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Therefore, 


(1  -  r\D)  n2  +  2n  +  1-fAjA.) 0  .  (118) 

Therefore, 

»  -  -3  *  (1  -  V  M/V  t  -g*  /  (1  -  y  ,  (119) 

from  which  <7 2/«fa  can  be  determined  from  the  calculations  already  available. 

From  Equation  (108), 

A2/Am  =  (VfnJ^)  (i  +  n2nD/{fru-l)2-n2r)  }  )*  (120) 

opt  u 

=  («  +  (AJAj)  n  2  }/  r)D}*  .  (121) 

0  opt 


Chapter  7 


COMPARISONS  BETWEEN  THEORY  AND  EXPERIMENT 
FOR  STATIC  EDUCTORS 

THE  LOCKHEED  EXPERIMENTS 


Reference  16  reports  a  comprehensive  test  program  on  two~dimer sional  eductors , 
similar  to  that  used  in  the  Lockheed  XV-4A  aircraft.  Because  of  t  he  large  size 
of  the  models  used,  the  measurements  can  be  anticipated  to  be  more  than  ade¬ 
quately  accurate. 

The  total  pressure  loss  through  the  duct  was  not  measured  separately,  so  that  it 
is  difficult  to  know  what  value  to  assign  to  the  loss  factor  np  .  However,  since 
both  the  diffuser  semi-angle  and  the  expansion  ratio  are  conservative.it  is  reason¬ 
able  to  expect  nr;  to  lie  in  the  range  0. 90  -  0. 95. 

The  arbitrary  geometry  theory  of  the  previous  chapter  has  been  worked  out  for  the 
geometry  of  the  Lockheed  eductors,  using  =  0.90  and  0.  95. 

As  shown  in  Figure  24,  the  agreement  obtained  is  good.  From  the  positions  of 
the  experimental  points  relative  to  the  theoretical  curves, the  actual  loss  factor 
r\D  can  be  deduced.  This  is  plotted  in  Figure  25.  This  variation  seems  rea¬ 
sonable. 

The  measured  entrained  air  ratio  (  n  )  is  seen  in  Figure  26  to  be  somewhat 
higher  than  the  theoretical  value,  a  result  which  may  be  attributable  to  the  non- 
uniform  flow  distribution  which  occurs  in  practice. 

In  general,  the  agreement  between  theory  and  experiment  appears  good,  and 
gives  confidence  that  "optimum  eductors"  can  be  built  with  much  better  efficien¬ 
cies.  On  the  Lockheed  eductors,  for  example,  a  simple  redesign  of  the  diffuser 
to  give  optimum  mixing  pressure  should  give  the  roughly  50  percent  increase  in 
augmentation  noted  in  Figure  27.  Improvement  in  diffuser  efficiency  w'ould  give  a 
further  increase,  and  should  not  be  hard  to  achieve  on  this  particular  configuration. 

A  SURVEY  OF  EXPERIMENTAL  RESULTS  GIVEN  IN  THE  LITERATURE 


Table  ID  lists  the  "best  results"  obtained  by  various  investigators.  That  is  to  say, 
when  an  investigator  experimentally  measured  the  variation  of  J g/d 'a  with  some 
parameter,  such  as  diffuser  area  ratio,  for  example,  Table  HI  lists  the  conditions 
under  which  he  obtained  the  maximum  value  of  J 2/Ja  •  ^  *8  usual  to  plot  results 
such  as  these  against  a  geometric  area  ratio,  rather  than  the  entrainment  ratio  ^  ), 
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DIFFUSER  AREA  RATIO  A^Am 


Figure  25.  Inferred  Variation  of  Effective  Diffuser  Efficiency 
With  Area  Ratio  for  the  Lockheed  Eductors. 
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Figure  26.  Predicted  and  Measured  Mass  Flow 
Ratio  n  for  the  Lockheed  Eductors. 


because  (  n  )  is  rarely  determined  during  test. 


For  parallel  tube  eductors,  Figure  28  shows  that  the  state  of  the  art  corresponds 
roughly  to 

J2/Ja  «  0.614  0.636  logJ0  (A/Aj)  ,  (122) 


where  A  is  the  eductor  tube  cross-sectional  area  (Aj  +  Afl  =  Ag  . 

The  highest  augmentation  was  reported  by  Payne  in  Reference  11.  It  was  achieved 
in  a  program  directed  by  N.  K.  Walker,  who  was  chiefly  responsible  for  the  excel¬ 
lence  of  the  results. 

It  is  of  interest  to  note  that  the  multi-staged  eductors  reported  by  Morrison  in 
Reference  2  appear  to  show  no  advantage  over  the  single  stage  when  the  eductor 
tubes  are  parallel-walled. 

When  the  duct  is  not  parallel,  the  appropriate  area  ratio  for  correlation  is 
Ag/A  j  .  As  shown  in  Figure  29,  the  "apparent  state  of  the  art"  is  now  some¬ 
what  better  for  smaller  entrainment  ratios ,  being  given  by 


=  0.995  +  0.565  log10  (Ag/Aj  .  (123) 

The  Payne  data  points  still  give  the  highest  augmentation  and  tend  to  depress  the 
apparent  state-of-the-art  line  below  its  true  position  at  the  higher  area  ratios,  be¬ 
cause  they  were  obtained  with  constant  area  tubes.  It  is  most  probable  that  educ¬ 
tors  fitted  with  diffusers  would  have  achieved  better  results  in  this  region  of  high 
area  ratio. 

In  Figure  30  the  data  points  are  compared  with  the  "optimum  eductor"  theory  of 
this  report.  They  correspond  to  equivalent  diffuser  efficiencies  in  the  range  0. 8 
<  <  0. 9.  Since  much  of  the  inefficiency  of  the  tested  eductors  must  be 

charged  to  non-optimum  mixing  pressure,  that  is,  incorrect  area  ratios,  their 
actual  internal  losses  will  be  less  than  Figure  30  implies.  Probably,  the  Lock¬ 
heed  eductor  losses  inferred  in  Figure  25  are  more  generally  representative. 

The  area  ratio  A^/A  •  is  sometimes  used  as  a  means  of  correlating  augmentation 
ratio.  As  shown  in  Figure  31,  this  does  not  reduce  the  scatter.  Since  it  is  a  less 
fundamental  parameter  than  A  2/Aj  ,  its  use  is  not  recommended. 


66 


AUGMENTATION  RATIO 


AREA  RATIO  A^/Aj 

Figure  27.  Performance  of  the  Lockheed  (Reference  1) 
Eductor  Geometry  as  a  Function  of  Diffuser 
Exit/Primary  Nozzle  Area  Ratio. 
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Figure  28.  Best  Augmentation  Ratios  Obtained  by  Various 
Investigators  with  Constant  Area  Eductors. 


APPARENT  STATE  OF  THE  ART 
FOR  MULTI-STAGE 
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Figure  29.  Best  Augmentation  Ratios  Obtained  by  Various 
Investigators,  as  a  Function  of  Exit  Area  Ratio. 


Optimum 


MULTI-STAGE 


Figure  31.  Best  Augmentation  Ratios  Obtained  by  Various 
Investigators,  as  a  Function  of  Mixing  Chambe 


Table  HI 

BEST  AUGMENTATION  RATIO  MEASURED  BY  VARIOUS  INVESTIGATORS 


SINGLE  STAGE  STATIC  EDUCTORS 


SYMBOL 

AUTHOR 

A2M2  Aj/Aj 

A2/Ay 

Wa 

REMARKS 

0 

R.  Morrison 

1.0 

15.0 

15.0 

1.35 

(Data  Abstracted 

(Reference  2) 

1.0 

20.0 

20.0 

1.37 

from  Reference 

1.0 

30.0 

30.0 

1.40 

17) 

1.0 

37.0 

37.0 

1.40 

0 

Peter  R.  Payne 

1.0 

130.5 

130.5 

1.95 

Axisymmetric. 

(Reference  11) 

1.0 

340.0 

340.0 

2.11 

1.0 

330.0 

330.0 

2.21 

No  Diffuser. 

1.0 

370.0 

370.0 

2.21 

• 

G.  L.  Rabeneck 

1.73 

52.0 

90.0 

1.89 

2-D  Models. 

P.  K.  Shunpert 

1.61 

30.3 

48.8 

1.83 

and  J.  F.  Sutton 

1.38 

13.0 

18.0 

1.50 

(Reference  16) 

1.23 

7.  55 

9.3 

1.36 

■ 

i  aan  Bertin 

_ 

— 

6.0 

1.20 

Presumably 

(Reference  17) 

— 

— 

10.0 

1.30 

Axisymmetric. 

— 

— 

14.0 

1.40 

— 

— 

20.0 

1.50 

© 

M.  F.  Gates  and 

2.0 

12.5 

25.0 

1.50 

2-D  Model  Test. 

J»  W.  Fairbanks 

1.0 

15.0 

15.0 

1.37 

0°  Diffuser  Angle. 

(Reference  18) 

— * 

— 

25.0 

1.49 

30°  Diffuser  Angle 

— 

— 

35.0 

1.245 

60°  Diffuser  Angle 

1.0 

22.5 

22.5 

1.42 

Axisymmetric. 

2.0 

11.8 

23.6 

1.56 

Axisymmetric. 

▲ 

K.  Cossairt 

3.11 

35.5 

111.0 

1.55 

Axi  symmetric. 

(Reference  19) 

2.73 

55.5 

153.0 

1.66 

Stalled  Diffuser 

(TCREC  62-66) 

2.36 

74.4 

177.0 

1.79 

Flow. 

1.36 

135.5 

189.0 

1.78 

A 

Peter  R.  Payne 

1.65 

14.4 

23.8 

1.39 

Axi  symmetric. 

and  Alastair 

2.02 

31.4 

63.5 

1.71 

Stalled  Diffuser 

Anthony 

2.27 

51.0 

115.8 

1.90 

Flow. 

(Reference  27) 
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♦ 

Chia-An  Wan 

1.0 

4.0 

4.0 

1.03 

Axi  symmetric. 

(Reference  29) 

1.0 

12.0 

12.0 

1.20 

1.0 

16, 0 

16.0 

1.26 

1.0 

36.0 

36.0 

1.40 

1.0 

47.0 

47.0 

1.48 

1.0 

105.0 

105.0 

1.50 

MULTI-STAGE  EDUCTORS 

SYMBOL 

AUTHOR 

No. 

of  Stages  A  t /A; 

J2/Jc 2 

REMARKS 

+ 

Reaves  Morrison 

3 

13.0 

1.40 

Data  Taken  From 

(Reference  2) 

f 

U 

31.0 

1.52 

Reference  17. 

3 

43.0 

1.56 

Constant  Section 

3 

70.0 

1.62 

Duct. 

X 

Jean  Bertin 

? 

35.0 

1.95 

Possibly  Three 

(Reference  17) 

? 

40.0 

2. 10 

Stages. 

? 

45.0 

2.00 

? 

50.0 

2.  30 

A2/Aj 

A 

Paul  Guienne 

? 

71.0 

1.92 

Possibly  Three 

(Reference  17) 

? 

84.0 

2.16 

Stages. 

? 

102.0 

2. 32 

Chapter  8 


SOME  SPECIAL  SOLUTIONS 

PERFORMANCE  OF  AN  EDUCTOR  IN  AN  AXIAL  STREAM 


Earlier  in  this  report  it  was  observed  that  the  augmentation  ratio  can  be  expect¬ 
ed  to  be  less  when  an  eductor  is  moving  axially  or  is  immersed  in  an  axial  stream. 
This  phenomenon  is  due  to  two  effects.  First,  the  secondary  airflow  entering  the 
eductor  corresponds  to  a  momentum  drag  of  nm^u0  ,  and  this  can  be  wholly  can¬ 
celled  only  by  additional  thrust  if  there  are  no  losses  in  the  eductor.  Since  loss 
must  occur,  and  the  energy  content  of  the  secondary  fluid  is  relatively  large,  the 
eductor  experiences  a  rapid  fall-off  in  thrust  with  increasing  speed. 

The  second  reason  is  more  fundamental.  Relative  to  vehicle  axes,  the  jet  has  a 
kinetic  energy  of  -u2: .  A  certain  ratio  (  )  of  this  can  be  transferred  to 

the  secondary  flow.  That  is . 

(uz2  -  u2o)  =  Ip  irhj-uj  ; 

therefore, 

( uJvl  -)2  =i'/n  +  (u  Ju  J2  .  (124) 

6  °  J 

But  the  thrust  gain  is  given  by 

*  -  n  [fa 2/"j)  ~  (uc/uj 


=  n[{6l>/n)  +  (u/u-)2)*-u  /u .  )  ,  (125) 

O  J  0  J 

which  tends  to (  -  n  (u/u-)  as  (u  /u-)  ■*  0 


Consequents,  as  the  axial  speed  increases,  the  same  energy  exchange  between 
the  primary  and  secondary  flows  is  progressively  of  less  value  in  producing  a 
thrust  increase. 
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Expressing  the  previously  given  general  eductor  theory  in  a  form  more  suitable 
for  studying  this  effect,  the  augmentation  ratio  is 


J?/Ja  *  {(*+Um.u2  -  nmju0)  /  (126) 

=  (”/Mja)2  (u2/uj}  (n  +  1}  -  (k/n3a)2  (u</u3}  n  (127) 

*  (VA^.}  [kp2  (n+l)2  +  T)D  (  1+n  {(uo/Uj)2-Lp2  )2 }%-  nuju-)  (128) 


Figure  32.  Augmentor  in  an  Axial  Stream. 


Taking  the  optimum  value  of  Ap^  given  by  Equation  (83)  , 

J 2/* 'o  ~  (f n+1 ) /AP  •)  ( { n^/ {(n+1)2 -n2  ^D^+(  vQ/ uj) 2  (129) 

opt  J 

~{n/(n+l )}( Uq/uj))  . 

From  Equations  (91)  and  (92) ,  therefore, 


(n+D/t&j  -  [(n+1  )2  -n2n^)2  (n+1)  /  ( { (n+1  )2  -n2  n^)2  -n2  n^)  .  (130) 

Thus,  Equations  (129)  and  (130)  give  an  explicit  variation  of  with 

the  velocity  ratio  u  /u  •  .  To  determine  variation  with  u  ,  u-  must  be 
calculated.  3 

Now  APj  =  A p2  +  %pu2.  ; 

therefore, 

h?u2-  =  A P;  [. 1 -  Apj/ APj  )  (131) 

=  AP^/APj  from  Equation  (132).  Therefore, 


W2  =  r&VAV  A?j 


(132) 


(uo/uJ2  =.[{(n+l)2~n2r]^}2-n2T\2  )/[(n+l  )2-n2r\j~,y 


(133) 


Thus,  Equations  (129),  (130), and  (133)  completely  determine  the  augmentation 
ratio  as  a  function  of  A P  / bPj  .  They  are  used  to  produce  the  curves  of  augmen¬ 
tation  ratio,  as  a  function  of  n  and  AP  /AP  -  in  Figures  33  -  36,  taking  a 
typical  loss  coefficient  of  =  0. 95. 

It  is  evident  that  a  significant  degree  of  augmentation  is  obtained  only  when  the 
total  pressure  ratio  AP^/APj  is  quite  low.  However,  there  may  be  applications 
where  this  requirement  is  met.  An  example  might  be  the  application  of  an  eductor 
to  water  jet  propulsion  of  a  ship  .where  the  eductor  intake  is  immersed  in  the  ship 
boundary  layer  so  that  its  effective  intake  velocity  is  much  lower  than  u0 
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OPTIMUM  AUGMENTATION  RATIO  'V*7, 


Figure  34.  Cross-Plot  of  Figure  33  for  Low 
Entrainment  Ratios. 


OPTIMUM  AUGMENTATION  RATIO 


0  .02  .04  .06  .08  .1 

TOTAL  PRESSURE  RATIO  APa/APj 


Figure  36.  Best  Possible  Augmentation  Ratio,  as  a 
Function  of  the  Pressure  Ratio  &(/ 

(  n d  =0. 95). 
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THE  TOTAL  HEAD  RISE  ATTAINABLE  THROUGH  AN  EDUCTOR 


Id  the  application  of  eductor  technology  to  many  problems,  the  total  head  rise 
through  the  eductor  may  be  of  greater  interest  than  its  thrust  or  force  augmen¬ 
tation.  The  Martin  recirculation  GEM  lift  system  (Reference  19)  is  a  particu¬ 
larly  good  example.  The  theory  of  earlier  chapters  of  this  report  is  accordingly 
reworked  in  this  section  to  include  an  analysis  of  total  pressure. 

Since  a pn  +  =  Sp.^  +  no  [l+n(Uj/uj)2/(n+l)2  ,  (134) 

the  total  pressure  rise  through  the  eductor  is 

tJi  *  AP^  -A P2  =  nD  hQU2-{{l+n(u1/uj)}2/(n+l)2)-  %p u2  (135) 

bH/hp  »?. 

=  {r\I/(n+l)2  y[2nx\D/(n+l)2)(u1/u^)  +  [rfinp/fn+l  )2-l](u1/iij)2  .  (136) 


Differentiating  with  respect  to  (  i^/iy  )  and  equating  to  zerj  to  find  the  value 
of  (  i f/uj  )  for  maximum  total  head  rise, 

2nr\jj/ (n+1  )2  =  2(1^/ Uj)[l-{n2T\p/(n+l  )2}) 


-■  n*\y/{  (n+1  )2-n2r\D) 


opt 


Substituting  Equation  (137)  in  Equation  (136), 


LE/%q 


*§|  *  ^((n+ljt-n2^) 

I  opt 


(137) 


(138) 
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rr 

4 

i; 


S 


* 


or 


(i/rcj  (nj/uj) 


opt 


(139) 


Note  that  since 


Apj  =  /%>  =  AP0-%pu|/%p  -u-)1 


bPQ  -  n2  n^/{  (n+1)2 -n2 }2 


(140) 


This  is  the  same  optimum  mixing  pressure  as  that  obtained  earlier  for  maximum 
thrust  augmentation.  Thus,  optimization  for  maximum  total  head  rise  gives  the 
same  geometry  as  optimization  for  maximum  thrust  augmentation. 

RECIRCULATING  EDUCTORS 


In  this  sectiou,  the  "optimum  eductor"  concept  is  extended  to  d  elude  the  effect 
of  recirculation.  Examples  of  recirculating  eductors  are  jet-driven  gas  or  water 
tunnels,  and  the  Martin  recirculating  GEM  system  (Reference  19). 

When  the  diffuser  exit  static  pressure  is  equal  to  ambient.  Equation  (84)  becomes 


J2/Ja 


opt 


the  second  term  being  the  effect 
head. 


Now, since 


+  (n+DtfbPc/bPJ  1/bPj  ,  (141) 

of  the  increase  above  ambient  of  the  inlet  total 


J2/Ja  (J2/Jo)  =i(n+2)/ APjO  (u2/uj) 

“2  =  (J2/Ja}  U/(n+1  j  > 


(142) 

(143) 


82 


OPTIMUM  TOTAL  HEAD  RISE  PARAMETER  nbH/hpu 


as* 


APg/AP^.  «  =  ^ 2/^ o}2  1/(Un)2 

&p  /ap.  =  €aPo/ap  .  =(  sap  ./(«+:; 2  )<v /«/■  ;2 

d  0  V  J  U  yl 

Substituting  in  Equation  (141), 


(144) 


(J2/Ja)2  =  (J2/Ja)2  +  t(J  /Ja)2  *  (145> 

therefore, 

=  (i/a-u%)rj  /j  )  (146) 

opt  ^  opt 

pc-°  * 

S  being  the  total  head  recovery  factor. 

Thus,  optimum  static  eductor  curves  can  be  used  directly  to  find  the  optimum 
augmentation  possible  with  recirculation. 

For  c  38  20%  40%  60%  80%  90% 

l/Ci-V*  *  1.117  1.29  1.58  2.235  3.16. 

In  other  words,  the  thrust  increase  due  to  recirculation  is  not  particularly  under 
optimum  conditions,  until  the  total  head  loss  falls  below  20  percent.  Such  a  low 
head  loss  is  difficult  to  achieve. 


s 
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Chapter  9 


SOME  EDUCTOR  COMPONENT  DESIGN  CONSIDERATIONS 


The  quality  of  the  eductor  design  has  so  far  been  considered  only  in  terms  of 
"effective  diffuser  efficiency"  (  np  ),  which  accounts  for  all  the  total  pressure 
losses  attributable  to  the  duct  or  shroud.  It  is  clear  that  the  value  of  (  nD  ) 
actually  achieved  will  depend  largely  upon  the  skill  with  which  the  duct  is  shaped. 
The  purpose  of  this  chapter  is  to  review  briefly  the  major  design  variables  which 
influence  overall  performance. 

It  should  be  recognized  that,  in  addition  to  mixing  losses,  there  are  four  main 
sources  of  energy  loss : 

1.  Intake  loss. 

2.  Loss  attributable  to  the  secondary  air  flowing  past  the  primary  nozzles. 

3.  Wall  skin-friction  loss  in  the  mixing  chamber. 

4.  Diffuser  loss. 

The  reference  velocity  for  the  first  three  of  these  can  conveniently  be  . 
Thus,  if  ( 7  -  n  rt  )  Ho  ^  is  the  loss  involved, 

{(l-r\j)  +  (1- r\2)  +  (1~ %PI^  *  (l~r\ DE)  Hpi£  ; 

therefore, 

tt-n De}  =  rW2  l  •  (147) 

From  Equation  (76)  , 

uj/\  =  (yl+l)/(l+n'ui/Uj-  ; 

therefore, 

Kj/%  =  ( =  (n+1  )/{(vj/iq}+n)  .  (148) 

This  is  plotted  in  Figure  38  for  a  typical  optimum  eductor.  Vfo**  that,  particu¬ 
larly  at  the  smaller  entrainment  ratios,  a  change  in  an  upst  ream  component  ef¬ 
ficiency  %  causes  a  significantly  smaller  change  in  the  effective  diffuser 
efficiency  n#  .  In  part,  then,  this  constitutes  some  justification  for  the 
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ENTRAINMENT  RATIO  n 


Figure  38.  Typical  Variation  of  the  Velocity  Ratio  (uj/ u^f) 
with  Entrainment  Ratio  for  a  Diffuser  Efficiency 
Op  =0.9. 
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simplification  introduced  by  referring  all  losses  to  the  diffuser  throat  position. 
THE  INTAKE  LOSS 


In  Reference  16,Rabeneck  et  al  considered  an  eductor  whose  bellmouth  intake 
lip  was  made  up  of  a  circular  arc.  After  testing  various  lip  sizes,  they  con¬ 
cluded: 

"An  entrance  radius  approximately  equal  to  the  ejector  mixing  section  width 
appears  to  be  optimum.  A  radius  of  0.40  times  the  mixing  section  width,  how¬ 
ever,  results  in  only  a  2. 5  percent  reduction  in  augmentation  factor. 

"Ejector  performance  is  sensitive  to  ejector  entrance  total  pressure  losses  in 
the  secondary  stream. " 

Other  experimenters  have  also  pointed  out  that  their  results  appear  sensitive  to 
both  the  size  and  the  shape  of  the  inlet,  although  without  suggesting  any  under¬ 
lying  reason.  The  following  arguments  are  presented  as  an  explanation  for  at 
least  some  of  the  observed  anomalies. 

When  a  fluid  flows  around  a  curved  surface,  the  centrifugal  force  field  generated 
in  the  flow  results  in  the  wall  static  pressure's  being  lower  than  it  would  be  in 
rectilinear  flow  under  the  same  conditions.  Among  other  references,  this  sub¬ 
ject  is  discussed  in  detail  by  Payne  and  Anthony  in  Reference  84. 

If  this  curved  surface  is  tangential  to  a  straight  surface,  the  fluid  near  the  wall 
must  experience  a  rapid  increase  in  static  pressure  when  it  passes  the  point  of 
tangency.  This  amounts  to  a  "rapid  diffusion"  from  a  velocity  Uj  to  a  lower 
velocity  Kg  ,  which  Payne  and  Anthony  have  shown  to  result  in  a  total  pressure 
loss. 


ts/ho  =  [i-( u2/i^))2 


(149) 


An  example  of  such  a  rapid  diffusion  just  behind  a  bellmouth  is  given  in  Figure 
39,  the  diffusion  taking  place  between  pressure  taps  8  and  10. 

It  follows  that  a  discontinuity  in  curvature  is  undesirable  in  an  inlet.  Also,  it 
i3  preferable  to  design  for  an  approximately  uniform  wall  pressure,  minimizing 
the  diffusion  necessary  at  the  bellmouth  exit, assuming  that  it  is  distributed  over 
a  finite  wall  length.  There  is  a  real  need  tor  both  theoretical  and  experimental 
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WALL  STATIC  PRESSURE  IN  LB.  /FT 
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- 1  . 


DISTANCE  ALONG  WALL  IN  INCHES 


Figure  39.  Static  Pressure  Along  an  Axiaymmetric 
Eductor  -  Payne,  Inc. ,  Data. 
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investigations  of  bellmouth  inlets  designed  to  such  criteria.  Almost  all  work 
to  date  has  been  ac  hoc  in  nature.  In  particular,  it  should  be  noted  that  although 
an  empirically  designed  bellmouth  may  have  small  losses  by  normal  standards, 
these  losses  may  still  cause  a  major  loss  of  thrust  augmentation  in  an  eductor 
application. 

The  smallest  bellmouth  is  likely  to  be  one  with  a  constant  wall  pressure.  Two- 
dimensional  solutions  exist  for  this  case:  the  free-boundary  Borda-mouthpiece 
solution,  and  the  free-boundary  solution  for  fluid  issuing  from  a  sharp-edged 
aperture.  From  Lamb  (Inference  85,  Art.  74),  the  solution  for  the  Borda- 
mouthpiece  is 


x  =(2b/v  )(sin2  %0  -  log  sea  %0  ) 

y  =(b/ nj(0  -  sin  Q)  , 


(150) 


where  2  b  =  final  distance  between  walls  . 

Also  from  Lamb,  Art.  75,  the  jet  issuing  from  a  sharp-edged  aperture  has  a 
boundary  defined  by 

x  =( 4b/n)  sin 2  %0 

y  =(2b/n)  {lag  tan  {(k/4)  +  (Q/2))  -  sin  0} 

0  <  0  <  %n 


The  first  of  these  solutions  is  appropriate  to  the  intake  of  a  static  eductor, 
where  the  entrained  air  can  flow  from  any  direction  between  0  and  180°  to  the 
eductor  axis.  The  second  solution  is  appropriate  for  an  intake  mounted  in  a 
plane  wall. 

Since  either  curve  is  asymptotic  only  to  a  parallel  duct,  it  might  seem  that  at 
some  point  the  theoretical  curve  must  be  empirically  faired  into  the  duct  wall  of 
the  mixing  chamber.  However,  in  order  to  maintain  constant  static  pressure,  it 
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is  necessary  to  progressively  reduce  the  area  of  the  mixing  chamber  moving 
downstream  from  the  nozzles.  This  compromise  can  therefore  be  minimized 
by  continuing  part  of  the  theoretical  intake  curve  into  the  mixing  chamber. 

The  free-streamline  solutions  given  earlier  in  this  section  are  for  two-dimen¬ 
sional  flow.  A  convenient  method  of  applying  them  to  the  axisymmetric  case 
is  to  assume  that  the  ordinate  y  is  proportional  to  the  square  of  the  radius. 
The  accuracy  of  this  method  is  not  known.  The  only  justification  for  its  use  is 
that  it  is  frequently  employed  and  often  gives  good  results. 

The  axisymmetric  Borda-mouthpiece  streamline  obtained  in  this  way  is  given 
in  Figure  40  and  Table  IV,  together  with  the  two-dimensional  solution  from  which 
it  was  derived. 

THE  PRIMARY  NOZZLE  LOSS 


The  secondary  air  flowing  past  the  primary  nozzles  suffers  a  total  pressure  loss 
which  corresponds  to  the  "drag"  force  which  it  develops  on  the  nozzles.  In  terms 
of  one-dimensional  flow,  the  total  pressure  loss  is  given  by 


4JAPI  =  (Ajfjjg)  .  (152) 

The  nozzle  drag  coefficient  will  be  higher  than  for  a  similar  body  in  an 

infinite  stream,  for  two  reasons.  In  the  first  place,  its  proximity  to  the  eductor 
walls  will  give  rise  to  a  "wind  channel  blockage"  effect.  Secondly,  the  "jet  drag" 
effect  identified  by  Payne  and  Anthony  in  Reference  30  will  give  rise  to  lower 
static  pressures  on  the  downstream  faces  of  the  primary  nozzle  array,  further 
increasing  its  apparent  drag  in  the  surrounding  secondary  flow.  The  magnitude 
of  both  these  effects  depends  upon  the  primary  nozzle  geometry. 

By  moving  the  primary  nozzles  out  from  the  throat  of  the  inlet  bellmouth,  they 
enter  a  region  of  lower  primary  flow  velocity  and  hence  reduce  the  magnitude  of 
these  losses.  However,  this  now  means  that  some  of  the  mixing  occurs  at  a 
greater  static  pressure  than  that  prevailing  in  the  mixing  tube,  so  that  the  educ¬ 
tor  efficiency  is  reduced.  Thus,  there  is  an  optimum  axial  position  for  a  given 
nozzle  array  and  bellmouth,  which  will  depend  upon  individual  geometrical  var¬ 
iables.  Any  simple  generalization,  such  as  some  investigators  have  been  tempted 
to  make  on  the  basis  of  isolated  tests,  is  not  necessarily  applicable  to  new  designs. 
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Figure  40.  The  Borda  Mouthpiece  Free  Streamline 
Solution  as  an  Intake  Profile. 


Table  IV 


THE  ’’BORDA -MOUTHPIECE"  SOLUTION 


e 

x/b 

y/b 

2  -y/b 

(2-y/b)h 

20 

0.0094 

0.0022 

1.9978 

1.4140 

40 

0.0349 

0.0175 

1. 9825 

1.4070 

60 

0.0676 

0.  0579 

1.9421 

1.3930 

80 

0.  0932 

0. 1312 

1.8688 

1.3670 

100 

0. 0921 

0.2420 

1.7580 

1.3250 

120 

0. 0362 

0. 3910 

1.6090 

1.2680 

140 

-0. 1210 

0.  5740 

1.4260 

1.1930 

150 

-0.2660 

0. 6750 

1.3250 

1.1500 

160 

-0.4970 

0.  7790 

1.2210 

1.1030 

165 

-0. 6700 

0. 8340 

1.1660 

1.0800 

170 

-0. 9200 

0.  8880 

1.1120 

1.0530 

175 

-1. 3600 

0.9430 

1.0570 

1.0285 

180 

_  CO 

1.0000 

1.0000 

1.0000 

THE  MIXING  CHAMBER  WALL  LOSS 


A  rough  indication  of  the  skin  friction  loss  in  the  mixing  chamber  can  be  ob¬ 
tained  by  assuming  that  the  secondary  inlet  velocity  is  applied  over  the  entire 
surface.  When  the  mixing  section  is  properly  designed  for  constant  static 
pressure,  this  provides  a  reasonably  accurate  estimate,  and  we  obtain 


P1  "  Pm  =  C  f  (SWET^\)  ’  ^153^ 

For  the  axisymmetric  case,  the  theory  developed  in  this  report  then  gives 

(Pj/P  V%pw&  =  C  f  (x/dm  )  {  ( n+l)2n 2  (AjMj  J2 >  (®  ( ^Aj/\  )+l  )*  (154) 

/  U+*2  A/A2  Y, 

where  x  =  the  mixing  chamber  length 
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=  diameter  of  the  diffuser  throat  . 

Equation  (154)  is  plotted  in  Figure  41.  For  a  typical  skin  friction  coefficient  of 
0. 003,  for  example,  and  a  mixing  length  ratio  of  x/'c^  =  5. 0,  Figure  41  gives 

an  equivalent  diffuser  loss  of  about  5  to 7  percent.  Thus,  the  effective  diffuser 
efficiency  would  only  be  about  93  to  95  percent,  even  if  the  diffuser  itself  were 
100  percent  efficient.  Therefore,  it  is  necessary  to  minimize  the  term  C  * 

(  x/4n  )  as  much  as  possible.  C  *  can  be  minimized  hy  lamlnarlzation  tenh- 
niques  designed  to  delay  turbulent  transition  of  the  boundary  layer  as  far  down¬ 
stream  as  possible.  C  *  may  be  minimized  by  using  a  multiplicity  of  primary 
nozzles.  A  circular  mixing  section  is  also  desirable  since  it  minimizes  the  wet¬ 
ted  area. 

THE  DIFFUSER  LOSS 


The  problem  of  minimizing  diffuser  loss  has  been  studied  exhaustively  by  many 
workers.  A  few  references  to  this  work  are  given  in  the  bibliography;  however, 
a  detailed  analysis  of  the  problem  would  be  inappropriate  in  this  report.  Figures 
42  and  43  illustrate  two  unconventional  ways  of  improving  diffuser  efficiency  in 
the  eductor  application.  They  may  be  of  value  in  indicating  that  the  almost  uni¬ 
versally  employed  conventional  configuration  for  eductors  may  not  be  optimum. 


ICTAL  PRESSURE  DROP  PARAMETER(YP2-PmV£  pu^/C^x/d^) 


ENTRAINED  AIR  RATIO  n 

Figure  41.  Effective  Diffuser  Loss  Due  to  Skin  Friction 
of  Mixing  Chamber  Wall  (Circular  Section). 
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Figure  42.  Impinging  Jet  Diffuser  With  Mixing. 

NOTE :  Impinging  jet  diffuser  has  constant  wall  pressure  during  the 
turn,  so  that  mixing  can  continue  efficiently  in  this  region,  so  long  as 
it  occurs  near  the  wall.  The  wall  pressures  are  favorable  to  boundary 
layer  stability  until  final  diffusion  has  commenced  outboard  of  station  X. 
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Figure  43.  Centrifugal  Diffuser  With  Mixing. 


NOTE :  Centrifugal  diffuser  is  designed  for  constant  wall  pressure, 
all  the  pressure  rise  occurring  at  the  exit.  Hence, efficient  mixing 
can  take  place  in  diffuser,  provided  it  is  close  to  the  wails. 
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Chapter  1C 


SOME  UNCONVENTIONAL  EDUCTORS 
MULTI-STAGE  EDUCTORS 


Figure  23  shows  that  eductor  efficiency  decreases  as  the  entrainment  ratio  (  n  ) 
increases  for  constant  effective  diffuser  efficiency.  Also,  the  optimum  diffuser 
area  ratio  increases  with  (  n  ),so  that  it  is  progressively  more  difficult  to  ob¬ 
tain  satisfactory  diffuser  efficiency  at  the  higher  entrainment  ratios  needed  to 
generate  worthwhile  augmentation  ratios. 

An  obvious  remedy  is  to  divide  an  eductor  into  several  "stages, "  the  total  efflux 
from  each  stage  constituting  the  "primary"  of  the  next,  and  larger,  stage.  This 
not  only  reduces  the  entrainment  ratio  per  stage,  but  also  reduces  the  diffuser 
area  ratio  requirement.  The  development  of  such  a  system  has  been  associated 
chiefly  with  Melot  and  Bertin  (Reference  17),  although  Morrison  (Reference  2) 
has  also  reported  some  multi-stage  experiments.  As  indicated  in  Table  HI  and 
Figure  29,  the  multi-stage  eductors  tested  have  developed  greater  augmentation 
than  the  best  single-stage  units,  as  indeed  should  be  expected  from  the  considera¬ 
tions  mentioned  above. 

It  has  been  shown  that  even  a  single-stage  eductor  is  far  too  complex  a  problem 
for  "cut  and  try"  engineering  to  be  effective  in  obtaining  high  efficiency.  Since 
multi-staging  increases  the  number  of  variables  by  an  order  of  magnitude,  it  is 
evident  that  the  design  of  optimum  staged  eductors  can  only  follow  extensive 
theoretical  analysis.  However,  the  potential  benefits,  in  terms  of  improved  per¬ 
formance,  are  quite  large.  It  underscores  a  requirement  for  initiation  of  appro¬ 
priate  research  work  in  this  area. 

COANDA  EDUCTORS 


"Coanda  flow"  or  "clingmg  flow"  can  occur  when  a  fluid  jet  issues  close  to  a 
solid  surface.  A  typical  example  is  illustrated  in  Figure  44. 

Since  the  outside  edge  of  the  jet  is  at  ambient  static  pressure  throughout,  the 
centrifugal  forces  generated  around  the  curve  must  be  balanced  by  a  lower  than 
ambient  static  pressure  at  the  curved  wall. 


AMBIENT 

PRESSURE 


- X  REGION  OF  RAPID  DIFFUSION 


SEPARATION  BUBBLE 


Figure  44.  A  Coanda  Flow. 

If  the  curve  intersects  with  a  tangential  straight  wall,  then  the  sudden  disappear¬ 
ance  of  the  centrifugal  forces  in  the  flow  must  result  in  a  sudden  static  pressure 
rise  at  the  intersection.  This  effect  is  amenable  to  the  rapid  diffusion  loss  analy  ¬ 
sis  mentioned  earlier  (Equation  149),  and  we  find  that  quite  large  pressure  losses 
are  to  be  expected  when  the  jet  is  thick.  Coupled  with  the  skin  friction  loss  and 
the  static  pressure  gradient  across  the  jet,  this  diffusion  loss  explains  why  lab¬ 
oratory  studies  of  the  Coanda  effect  do  not  show  any  thrust  augmentation  effect, 
even  though  the  inner  portions  of  the  jet  entrains  ambient  air  at  a  lower  than 
ambient  static  pressure. 

It  has  been  shown  that  thrust  augmentation  occurs  when  a  jet  mixes  with  ambient 
air  at  a  static  pressure  which  is  lower  than  ambient,  and  that  there  is  an  optimum 
static  pressure  for  maximum  augmentation;  mixing  at  pressures  away  from  this 
optimum  results  in  less  than  the  maximum  possible  augmentation. 
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Figure  45.  Variation  of  Augmentation  Ratio  J 2/Ja 

With  the  Mixing  Pressure  Parameter  Ap,  . 


In  Coanda  flow,  the  mixing  pressure  varies  from  ambient  (a p2  =  0,  point  A  in 
Figure  45)  at  the  outside  of  the  jet  to 


APjj  =(-t/ro)(24t/ro  ))/  (1+  t/rQ  }2  (155) 

( WALL 

at  the  wall  (point  B  in  Figure  45).  Thus,  the  mixing  pressure  can  be  optimum 
at  only  one  streamline,  and  thus  the  overall  augmentation  will  be  less  than 
optimum. 

The  second  factor  determining  augmentation  is  the  efficiency  at  which  the  flow  is 
diffused  back  to  ambient  from  the  mixing  pressure  Ap^  . 

It  is  obvious  from  Figure  18  that  worthwhile  thrust  increases  can  be  obtained 
only  when  the  diffuser  efficiency  is  high,  90  percent  or  greater.  For  rel¬ 
atively  thick  jets,  the  diffusion  loss  at  the  end  of  a  Coanda  curve  can  be  quite 
large,  with  as  low  as  0, 5,  so  that  high  augmentation  ratios  can  be  obtained 
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JET  THICKNESS  PARAMETER  t/rfJ 

Figure  46.  Variation  of  Entrainment  Ratio  with  Jet 

Thickness  for  a  Right- Angle  Coanda  Bend 
{Based  on  Equation  58). 
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only  by  keeping  the  jet  thin.  However,  the  losses  due  to  skin  friction  then  be¬ 
come  important. 

The  amount  of  ambient  fluid  entrained  in  the  jet  is  another  reason  for  using 
thin  jets  when  augmentation  is  required.  The  attainable  augmentation  ratio 
naturally  increases  with  the  amount  of  fluid  entrained  and,  as  indicated  in 
Figure  46,  this  becomes  important  only  below  t/vQ  =  0. 1.  Very  high  entrain¬ 
ment  ratios,  therefore,  require  a  jet  thickness  which  is  only  a  few  percent  of 
the  radius  of  curvature;  and  in  this  case,  boundary  layer  effects  start  to  be¬ 
come  important. 

Although  nothing  can  be  done  about  Coanda  mixing  at  nonoptimum  pressures,  we 
can  avoid  the  efficiency  loss  due  to  diffusion,  as  indicated  in  Figure  47. 

OTHER  NOVEL  EDUCTORS 


Since  a  booster  rocket  operates  mainly  within  the  Earth's  atmosphere,  and  is 
payload-limited  by  its  static  thrust,  the  possibility  of  air  augmentation  naturally 
arises.  Work  is  proceeding  on  both  static  eductors  for  low  speeds  and  "after¬ 
burning"  with  air  at  higher  speeds,  as  indicated  by  Avery  and  Dugger  in  Refer¬ 
ence  26. 

Although  this  report  is  concerned  with  the  traditional  eductor  concept  involving 
steady-state  viscous  mixing  between  primary  and  secondary  flews,  it  would  be 
appropriate  to  notice  a  number  of  interesting  new  closely  related  developments 
in  non-steady  flow. 

In  addition  to  viscous  shearing,  a  jet  can  be  used  to  mechanically  move  the 
secondary  fluid.  By  use  of  a  pulsating  jet,  Lockwood  and  Patterson  (References 
14  and  22)  have  measured  substantially  higher  augmentation  ratios  than  those 
in  steady-state  operation.  Conceptually,  they  generate  a  succession  of  "air 
pistons"  which  augments  the  effect  of  viscous  shear  in  "driving"  the  secondary 
fluid  along  the  eductor  tube. 

Foa  (Reference  32)  achieved  an  analogous  effect  by  arranging  for  primary  Jets 
to  issue  from  a  rotating  body,  so  that  they  describe  a  helical  path,  as  illustrated 
in  Figure  48.  Thus,  apart  from  the  viscous  mixing  effects,  the  jet  sheets  tend 
to  "mechanically"  move  the  air  as  would  a  fan.  In  the  words  of  Foa, 

"Crypto-steady  pressure  exchange  is  a  mode  of  direct  energy  transfer  between 
flows,  based  on  the  principle  that  two  adjacent  streams  which  are  both  isoener- 
getic  hi  the  same  frame  of  reference  will,  in  general,  exchange  mechanical 
energy  in  any  other  frame.  The  efficiency  of  this  process  is  potentially  high, 
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ENTRAINED  AIR 


Figure  47.  A  Coanda  Thrust  Augmentor. 

NOTE:  Due  tc  the  static  pressure  variation  across  the  jet,  the  mixing 
pressure  can  be  optimum  at  only  one  streamline.  To  offset  this,  diffusion 
losses  can  be  half  those  for  an  equivalent  (two-wall)  diffuser.  So  far  as  is 
known,  no  Coanda  devices  so  far  tried  has  used  a  diffusion  section,  so  that 
from  eductor  theory,  we  should  not  expect  augmentation  to  be  obtained. 
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because  a  change  of  frame  of  reference  is  reversible,  and  the  associated  trans¬ 
fer  of  energy  is  therefore  nondissipative. 

fAm  application  of  this  principle  to  thrust  or  lift  generation  is  discussed  for  the 
purpose  of  illustration,  Li  this  application  the  interacting  flows  are  steady  and 
isoenergetic  in  a  rotating  frame  of  reference  but  exchange  energy  in  a  stationary 
frame.  The  exchange  mechanism  is  essentially  similar  to  that  of  a  turbofan,  but 
the  ^blades  *  are  now  patterns  rather  than  bodies  of  abiding  material. " 

A  system  somewhere  between  those  of  Lockwood  and  Foa  is  the  oscillating  jet  of 
Saunders  (Reference  33).  Here  the  "mechanical"  interaction  between  the  jet  and 
the  secondary  fluid  is  similar  to  the  propulsion  obtained  by  a  fish  as  it  swims. 

By  using  liquid  amplifier  techniques  to  oscillate  the  primary  jet,  Saunders  avoids 
the  need  for  the  moving  mechanical  parts  required  by  the  Foa  system. 
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CONCLUSIONS 


This  report  reviews  the  state  of  the  art,  so  far  as  steady  flow  eductor  technol¬ 
ogy  is  concerned,  with  particular  reference  to  eductors  which  augment  the  thrust 
of  a  primary  jet.  Since  the  steady  flow  jet  pump  principle  can  be  applied  in  a 
number  of  ways,  and  since  non-steady  flow  systems  are  also  of  interest,  the 
total  field  is  much  wider  than  that  covered  by  this  review- 

So  far  as  steady  flow  eductors  are  concerned,  the  following  conclusions  are 
drawn. 


It  seems  clear  from  the  review  of  the  literature  that  a  great  deal  of  research 
work  remains  to  be  done  before  either  adequate  predictive  ability  is  achieved 
or  there  is  confidence  in  obtaining  something  near  maximum  efficiency. 

It  seems  equally  clear  that  hardware  applications  of  the  eductor  principle  have 
not  profited  to  the  extent  made  possible  by  existing  research  work,  despite  its 
limitations.  Failure  to  profit  from  Helmbold's  demonstration  of  the  advantages 
of  constant  mixing  pressure  is  a  case  in  point. 


The  existing  state  of  the  art  is  given  by 


Total  Augmentor  Thrust 
Primary  Thrust 


0.995  +  0.565  leg 2  Q  (A^AJ  , 


where  A2  and 


are  the  eductor  exit  and  primary  jet  areas  respectively. 


In  the  present  report, a  method  of  calculating  eductor  performance  is  presented 
and  is  shown  to  give  good  agreement  with  experiment.  Significant  improvements 
in  eductor  performance  can  therefore  be  expected  from  the  application  of  this 
theory  to  the  optimization  of  the  new  designs. 


Multi-stage  eductors  have  demonstrated  significantly  better  performance  loads 
than  equivalent  single-stage  units.  The  existing  state  of  the  art  is  difficult  to 
define,  however,  because  of  the  relative  porosity  of  experimental  data.  No  theo¬ 
retical  analysis  exists  for  the  case  of  multi-stage  units,  so  far  as  is  known,  al¬ 
though  there  would  be  no  difficulty  in  extending  the  theory  of  the  present  report 
to  this  more  complicated  case. 
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